Ejemplo n.º 1
0
Archivo: xxx.c Proyecto: Gilles86/afni
void elvis_CB( Widget w , XtPointer cd , XtPointer cb )
{
   XmDrawingAreaCallbackStruct * cbs = (XmDrawingAreaCallbackStruct *) cb ;
   XExposeEvent * ev = (XExposeEvent *) cbs->event ;
   int start , end ;
   static int ncall = 0 ;

   if( cbs->reason != XmCR_EXPOSE || ev->count > 0 ) return ;

STATUS("enter elvis_CB") ;

   if( need_plot ){
      make_plot() ;
      set_X11_background( XtDisplay(w) , XtWindow(w) , 255,255,255 ) ;
   }
   XClearWindow( XtDisplay(w) , XtWindow(w) ) ;

#if 1
   if( ncall%2 == 0 ) set_active_memplot( "Elvis" ) ;
   else               set_active_memplot( "Pelvis" ) ;
   ncall++ ;
#endif

STATUS("XClearWindow") ;

   start = end = 0 ;
#if 0
   memplot_to_X11( XtDisplay(w) , XtWindow(w) ) ;
#else
   memplot_to_X11_free( XtDisplay(w) , XtWindow(w) ) ;
#endif
   return ;
}
Ejemplo n.º 2
0
int main(int argc, char *argv[])
{
RGBPic *p;
PLOT *plot;

p=make_RGBPic(800,600);
RGBPic_clear_area(p,COLOR(60,80,50),0,0,p->width,p->height);

plot=make_plot(640,480);
plot->lower_x=3.1415;
plot->upper_x=12423;
plot->lower_y=2.34E-10;
plot->upper_y=43.5E-8;
draw_grid(p,plot,50,50);
free_plot(plot);

//RGBPic_printf(p,30,30,COLOR(255,240,210),COLOR(1,2,3),"Hello all !\n line 2 %d %p",argc,p);
//RGBPic_printf_v(p,300,300,COLOR(255,240,210),COLOR(1,2,3),"Hello all 2 !\n line 2 %d %p",argc,p);
RGBPic_draw_line(p,COLOR(40,255,80), 40,50,60,70);
RGBPic_draw_line(p,COLOR(250,255,80), 30,50,60,30);

RGBPic_dump_ppm("test1.ppm",p);
free_RGBPic(p);
return 0;
}
Ejemplo n.º 3
0
int main(int argc, char** argv)
{
	// --------------------------------------------------------------------
	// get job files and program options
	std::vector<std::string> vecScans;
	std::string strOutFile;
	bool bFlip = 0;
	std::string strVec0, strVec1;
	t_vec vec0, vec1;

	// normal args
	opts::options_description args("Extracts scan positions into a figure");
	args.add(boost::shared_ptr<opts::option_description>(
		new opts::option_description("scan-file",
		opts::value<decltype(vecScans)>(&vecScans), "scan file")));
	args.add(boost::shared_ptr<opts::option_description>(
		new opts::option_description("flip-y",
		opts::bool_switch(&bFlip), "flip y axis")));
	args.add(boost::shared_ptr<opts::option_description>(
		new opts::option_description("out-file",
		opts::value<decltype(strOutFile)>(&strOutFile), "output file")));
	args.add(boost::shared_ptr<opts::option_description>(
		new opts::option_description("vec0",
		opts::value<decltype(strVec0)>(&strVec0), "first scattering plane vector")));
	args.add(boost::shared_ptr<opts::option_description>(
		new opts::option_description("vec1",
		opts::value<decltype(strVec1)>(&strVec1), "second scattering plane vector")));

	// positional args
	opts::positional_options_description args_pos;
	args_pos.add("scan-file", -1);

	opts::basic_command_line_parser<char> clparser(argc, argv);
	clparser.options(args);
	clparser.positional(args_pos);
	opts::basic_parsed_options<char> parsedopts = clparser.run();

	opts::variables_map opts_map;
	opts::store(parsedopts, opts_map);
	opts::notify(opts_map);

	if(argc <= 1)
	{
		std::ostringstream ostrHelp;
		ostrHelp << "Usage: " << argv[0] << " [options] <scan-file 1> <scan-file 2> ... <scan-file n>\n";
		ostrHelp << args;
		tl::log_info(ostrHelp.str());
		return -1;
	}

	if(vecScans.size() == 0)
	{
		tl::log_err("No scan files given.");
		return -1;
	}


	if(strVec0.size())
		vec0 = tl::str_to_var<decltype(vec0)>("[3](" + strVec0 + ")");
	if(strVec1.size())
		vec1 = tl::str_to_var<decltype(vec1)>("[3](" + strVec1 + ")");
	// --------------------------------------------------------------------


	// first scan file serves as reference
	std::unique_ptr<tl::FileInstrBase<t_real>> ptrScan(
		tl::FileInstrBase<t_real>::LoadInstr(vecScans[0].c_str()));
	if(!ptrScan)
	{
		tl::log_err("Invalid scan file: \"", vecScans[0], "\".");
		return -1;
	}


	// get scattering plane if not already given as program arg
	if(vec0.size() != 3)
	{
		auto arrVec0 = ptrScan->GetScatterPlane0();
		vec0 = tl::make_vec<t_vec>({arrVec0[0],arrVec0[1],arrVec0[2]});
	}
	if(vec1.size() != 3)
	{
		auto arrVec1 = ptrScan->GetScatterPlane1();
		vec1 = tl::make_vec<t_vec>({arrVec1[0],arrVec1[1],arrVec1[2]});
	}

	if(bFlip) vec1 = -vec1;
	tl::log_info("Scattering plane: ", vec0, ", ", vec1, ".");


	std::vector<t_vec> vecAllHKL, vecAllPos;
	vecAllHKL.reserve(argc-2);
	vecAllPos.reserve(argc-2);

	// first scan position
	auto vecPos = get_coord<t_vec, t_real>(vec0, vec1, *ptrScan);
	tl::log_info(vecPos.first, " rlu  ->  ", vecPos.second, ".");
	if(vecPos.second.size() != 2)
	{
		tl::log_err("Invalid scan position for file \"", vecScans[0], "\".");
		return -1;
	}
	vecAllHKL.push_back(vecPos.first);
	vecAllPos.push_back(vecPos.second);


	// load other files
	for(int iFile=1; iFile<vecScans.size(); ++iFile)
	{
		std::unique_ptr<tl::FileInstrBase<t_real>> ptrScanOther(
			tl::FileInstrBase<t_real>::LoadInstr(vecScans[iFile].c_str()));
		if(!ptrScanOther)
		{
			tl::log_err("Invalid scan file: \"", vecScans[iFile], "\".");
			continue;
		}

		// remaining scan positions
		auto vecPosOther = get_coord<t_vec, t_real>(vec0, vec1, *ptrScanOther);
		tl::log_info(vecPos.first, " rlu  ->  ", vecPos.second, ".");
		if(vecPosOther.second.size() != 2)
		{
			tl::log_err("Invalid scan position for file \"", vecScans[iFile], "\".");
			continue;
		}

		vecAllHKL.push_back(vecPosOther.first);
		vecAllPos.push_back(vecPosOther.second);
	}


	if(!make_plot(strOutFile, vec0, vec1, vecAllHKL, vecAllPos, bFlip))
	{
		tl::log_err("No scan positions found.");
		return -1;
	}

	return 0;
}