comp.lang.idl-pvwave archive: archive » Help with code for hyperspectral image analysis

Home » Public Forums » archive » Help with code for hyperspectral image analysis

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: Help with code for hyperspectral image analysis [message #80528 is a reply to message #80374]

Wed, 20 June 2012 13:24

Kel
Messages: 7
Registered: June 2012

Junior Member

Update: Still running out of memory. I've already changed the values in the idlde.ini file, as recommended in other posts, but now I'm stuck. Google is suggesting that I increase the virtual memory of my pc, so I will give that a try.

On Wednesday, June 20, 2012 3:59:36 PM UTC-4, Kel wrote:
> Thanks you both so much! I had just gotten an out-of-memory error on Klemen's code, which seemed to be going ok up until that point, so I will try this modification now and hope it helps.
>
> On Wednesday, June 20, 2012 3:55:08 PM UTC-4, Brian J. Daniel wrote:
>> See my adjustments to Klemen's code below. This way you don't have to load the entire image, masked or otherwise, into memory. Also, untested.
>>
>> pro envi_masking_image
>>
>> ;select your image file
>>
>> ;Get your file
>> image=dialog_pickfile()
>>
>> ; Pick an output filename
>> output = dialog_pickfile()
>>
>> ;check to see if valid for ENVI
>> envi_open_file, image, r_fid=fid
>> IF fid[0] EQ -1 THEN Message,'Invalid file name: '+image
>>
>> ;learn about the image
>> ENVI_FILE_QUERY, fid, ns=ns, nl=nl, nb=nb, data_type=data_type,Interleave=interleave
>> help, data_type
>>
>> ;build array based on data type
>> ;data_holder=Make_Array(ns, nb, nl, Type=data_type) ; memory intensive
>>
>> ;make a mask
>> f = ENVI_GET_DATA(fid=fid, dims=dims, pos=140)
>> indx = where(f LT 40, count)
>> f[indx] = 0 ; now f is your mask
>>
>> ;apply the mask to all bands and save
>> OpenW, lun, output, /Get_Lun
>> if count gt 0 then begin
>> for i=0,nl-1 do begin
>> data = ENVI_GET_DATA(fid=fid, dims=dims, pos=i)
>> writeu,lun,data*f
>>
>> endfor
>> endif
>>
>> ; close binary file
>> close,lun
>> free_lun, lun
>>
>> ; write out an envi header file
>> envi_setup_head, Data_Type=data_type, /Open, R_FID=masked_fid, $
>> /WRITE, FNAME=output, Interleave=2 ; BIP
>>
>> end

Report message to a moderator

[Message index]

		Help with code for hyperspectral image analysis By: Kel on Mon, 11 June 2012 07:59
		Re: Help with code for hyperspectral image analysis By: Kel on Fri, 22 June 2012 06:51
		Re: Help with code for hyperspectral image analysis By: Brian Daniel on Fri, 22 June 2012 04:59
		Re: Help with code for hyperspectral image analysis By: Klemen on Fri, 22 June 2012 00:44
		Re: Help with code for hyperspectral image analysis By: Kel on Thu, 21 June 2012 19:42
		Re: Help with code for hyperspectral image analysis By: Klemen on Thu, 21 June 2012 14:44
		Re: Help with code for hyperspectral image analysis By: Kel on Thu, 21 June 2012 13:37
		Re: Help with code for hyperspectral image analysis By: Kel on Wed, 20 June 2012 13:24
		Re: Help with code for hyperspectral image analysis By: Kel on Wed, 20 June 2012 12:59
		Re: Help with code for hyperspectral image analysis By: Brian Daniel on Wed, 20 June 2012 12:55
		Re: Help with code for hyperspectral image analysis By: Klemen on Wed, 20 June 2012 10:57

Previous Topic:	Re: read irregular data
Next Topic:	sun symbol

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Mon Dec 01 02:48:46 PST 2025

Total time taken to generate the page: 3.19814 seconds