bouhennache rafik wrote: I try to import the ROI classes to my scatterplot but it's not easy my script IDL is ; Read in a 2-band image. file = FILE_WHICH('rad_2010.tif') !null = QUERY_IMAGE(file, info) gas = READ_IMAGE(file) red_band_pixels = REFORM(gas[4,*,*], PRODUCT(info.dimensions)) green_band_pixels = REFORM(gas[5,*,*], PRODUCT(info.dimensions)) ; Plot red versus green band. myPlot = SCATTERPLOT(red_band_pixels, green_band_pixels, $ SYM_FILLED = 5, SYM_COLOR = 'red', $ XTITLE = 'Red band value', $ YTITLE = 'Green band value', $ TITLE = 'Red vs. green band of "glowing_gas.jpg"') x = red_band_pixels y = green_band_pixels myLine = REGRESS(x, y, SIGMA=sigma, CONST=const, $ MEASURE_ERRORS=measure_errors) PRINT, myLine PRINT, 'Coefficients: ', myLine[*] end help

Could you clarify your question? It isn't clear what you are having trouble with or any error messages you are seeing.

Thank's Mr Zachary Norman for your reply. I have made a scatterplot and I want to import the four classes classification: urban soil vegetation watter to imporve the classes separation like this: I have the vector ROIs but it is difficult to me to write a script IDL to import it. your's