pro read_nadircam, infile set_plot, 'X' device,true_color=24,decompose=0,retain=2 loadct,39,/silent !P.Background=255 hdfid = HDF_SD_Start(infile, /read) hdf_sd_getdata,hdf_sd_select(hdfid,hdf_sd_nametoindex(hdfid,'camera_RGB')),cam_rgb hdf_sd_getdata,hdf_sd_select(hdfid,hdf_sd_nametoindex(hdfid,'camera_latitude')),cam_lat hdf_sd_getdata,hdf_sd_select(hdfid,hdf_sd_nametoindex(hdfid,'camera_longitude')),cam_lon hdf_sd_getdata,hdf_sd_select(hdfid,hdf_sd_nametoindex(hdfid,'DC8nav_latitude')),dc8_lat hdf_sd_getdata,hdf_sd_select(hdfid,hdf_sd_nametoindex(hdfid,'DC8nav_longitude')),dc8_lon hdf_sd_getdata,hdf_sd_select(hdfid,hdf_sd_nametoindex(hdfid,'DC8nav_ir_ground_temperature')),dc8_irtemp HDF_SD_End,hdfid window, 0,xsize=1600,ysize=900 region=[min(cam_lat), min(cam_lon), max(cam_lat), max(cam_lon)] if((region[2]-region[0])/(region[3]-region[1]) lt 9./16.) then begin mlat = 0.5*(region[0]+region[2]) region[0] = mlat-4.5/16*(region[3]-region[1]) region[2] = mlat+4.5/16*(region[3]-region[1]) endif else begin if((region[2]-region[0])/(region[3]-region[1]) gt 9./16.) then begin mlon = 0.5*(region[3]+region[1]) region[1] = mlon-8/9.*(region[2]-region[0]) region[3] = mlon+8/9.*(region[2]-region[0]) endif endelse map_set, 0,0,limit=region,/cylindrical,/isotropic,color=0,position=[0.02,0.02,0.98,0.98] nscans = size(cam_lat) nscans=nscans[1] for i=0l,nscans-1 do begin for j=0,1279 do begin tvlct,cam_rgb[0,i,j], cam_rgb[1,i,j], cam_rgb[2,i,j],0 oplot, [cam_lon[i,j]], [cam_lat[i,j]], psym=3,color=0 endfor endfor tvlct,0,0,0,0 for i=0l,n_elements(dc8_lat)-2 do oplot, dc8_lon[i:i+1], dc8_lat[i:i+1], color=255*(dc8_irtemp[i]-250.)/(280.-250.),thick=4 stop end