PRO RegisterGeoStationaryImage, filename ; File 200505030600x11_g.jpg obtained from: ; ; http://www.sat.dundee.ac.uk/pdus.html ; Need a filename? IF N_Elements(filename) EQ 0 THEN filename = Dialog_Pickfile(Filter='*.jpg') IF filename EQ "" THEN RETURN ; Read the image. Read_JPEG, filename, image Window, 2, XSize=625, YSize=625, Title='Original Image' TV, Rebin(image, 625, 625) PlotS, [12, 12, 613, 613, 12], [15, 611, 611, 15, 15], /Device, Color=Fsc_Color('red') ; Trim the image to just the image parts. (Found coordinates by inspection.) image = image[24:1227, 30:1219] ; Get size of image, create window. s = Size(image, /Dimensions) Window, 1, XSize=s[0]/2, YSize=s[1]/2, Title='Trimmed Image with Map Coordinates' ; Display image. LoadCT, 0 TV, Rebin(image, s[0]/2, s[1]/2) ; Set up map projection. (Found by inspection.) Map_Set, /Satellite, SAT_P=[6.5,0.0,0.0], XMargin=0, $ YMargin=0, /NoErase, 0.0, 63, $ Limit=[0.0, -5, 68.0, 63, 0.0, 131, -68.0, 63] ; Test by mapping continents. Map_Continents, Color=Fsc_Color('red') END