|
|
|
|
|
|
0 Replies and 3038 Views
Making a spatially resolved BPT 3038 0
Started by Deleted User
Hello, I'm trying to make a spatially resolved BPT. The problem is that i need to make an 'if' parameter to define a specific region. Something like that: mask=mask*0 for i=0, size(f&91;*,0,0&93;,/n_elements)-1, 1 do begin for j=0, size(f&91;0,*,0&93;,/n_elements)-1, 1 do begin if alog10(OIII&91;i,j&93;/Hb&91;i,j&93;) gt 0.0 then begin mask&91;i,j&93;=1 endif endfor endfor Those OIII and Hb are the y axis. The x axis is NII/Ha, so this ...
|
|
|
|
0 |
3038 |
19 May 2017 01:23 PM |
|
1 Replies and 2231 Views
Hypercube correction to measured spectrum 2231 1
Started by Deleted User
Hi, I am running ENVI 5.3 IDL 8.5 (64bit) I have hyperspectral images (.bil) of 9x9cm panels and each image also includes a Spectralon 50 reflectance standard. I want to convert the entire cube from DN to reflectance by scaling each pixel to the mean spectrum of the 50 reflectance standard. Any advice as to the easiest way to do this I right now I am thinking - select an ROI for the reflectance standard in my image, quick stats to get mean spectrum for the ROI, export mean spectrum to ...
|
|
|
|
1 |
2231 |
by MariM 16 May 2017 11:37 AM |
|
|
3 |
5491 |
by MariM15 May 2017 02:09 PM |
|
0 Replies and 3206 Views
Question about IDL Legend 3206 0
Started by Deleted User
Question from customer: Is there any way to mix a line and a single symbol in a legend You can reduce the number of symbol by setting the SAMPLE_WIDTH to zero, but then the line will disappear. Possible Answer: Unfortunately, I think the only way to reduce the number of symbols is by adjusting the SAMPLE_WIDTH. You can then use the POLYLINE routine to manually draw the line onto the legend: IDL> theory = SIN(2.0*FINDGEN(201)*!PI/25.0)*EXP(-0.02*FINDGEN(201)) ...
|
|
|
|
0 |
3206 |
25 Apr 2017 03:10 PM |
|
0 Replies and 2582 Views
Pressing down button doesn&39;t generate events 2582 0
Started by Deleted User
Question from user:
I am displaying a group of buttons with CW_BGROUP. I print out events when they happen. I don't see any event when I press the button but I see one when I release it. I do not have the /norelease
argument because I want to see the release event as well. I would like to see the event that occurs (if it occurs) when the button is pressed.
Possible answer:
Use the PUSHBUTTON_EVENTS keyword. An example is shown below:
pro dj_cas_19437_l1k4p5_event, event
COMPILE_OPT idl...
|
|
|
|
0 |
2582 |
06 Apr 2017 03:34 PM |
|
1 Replies and 3153 Views
Polar 3D 3153 1
Started by Deleted User
Can IDL do a 3d plot in polar coordinates Or do I need to convert to cartesian coordinates first
|
|
|
|
1 |
3153 |
06 Apr 2017 01:23 PM |
|
1 Replies and 3685 Views
Strange error with idlpy 3685 1
Started by Deleted User
Hi all,
here the strange error. I premise: two days ago i had no problem with my code.
Traceback (most recent call last):
File "G:/SVN EO/Simona/illegal_building_detection/Versione2/core/ibd.py", line 5, in
from idlcodes import *
File "G:\SVN EO\Simona\illegal_building_detection\Versione2\core\idlcodes.py", line 1, in
from idlpy import IDL
File "C:\Program Files\Harris\ENVI54\IDL86\lib\bridges\idlpy.py", line 318, in
class IDL(IDLTyp...
|
|
|
|
1 |
3685 |
03 Apr 2017 09:27 AM |
|
1 Replies and 3405 Views
Write_tiff multiband image 3405 1
Started by Deleted User
Hey all,
I am a student learning IDL and am seeking help in creating a multiband image from 11 .tiff files. I am trying to create a time sequence NDVI file that has each different NDVI time as a separate layer. I am having trouble building a code in IDL to make this happen.
I've been through the IDL help page for Write_tiff and have not been able to successfully construct a set of code to do what I want.
So I need to pull 11 images from a folder and combine them into one mu...
|
|
|
|
1 |
3405 |
24 Mar 2017 11:30 AM |
|
0 Replies and 3079 Views
idlpy management 3079 0
Started by Deleted User
Hi all,
i have a python script that calls and executes some python function making idl processes following the IDL.run('instruction') syntax.
More in detail, i have a set of function making "radiometric correction", "pansharpening" "layer stacking", and so on.
I make the import of idlpy module ad the begin of my file. The problem is that license remains busy throughout the execution duration.
I try also to make the import only from the function.
There is a way to im...
|
|
|
|
0 |
3079 |
14 Mar 2017 07:52 AM |
|
2 Replies and 2623 Views
Can&39;t standardize COLORBAR across multiple VECTOR maps 2623 2
Started by Deleted User
I'm creating wind maps using the MAP, VECTOR, and COLORBAR and LEGEND functions, and I'm using them to build wind animations. THE PROBLEM: I cannot standardize the color bar and sample vector length across multiple wind maps.
EXAMPLE-MAP-1
EXAMPLE-MAP-2
I'm using the TARGET option with COLORBAR to set the range:
cbar = COLORBAR(TARGET=wind_vectors, blah)
... so, I cannot use the RANGE=&91min,max&93 option with COLORBAR.
Also, the LEGEND sample vector change size from map to map. I'm u...
|
|
|
|
2 |
2623 |
11 Mar 2017 06:53 AM |
|
1 Replies and 2314 Views
Raster subset 2314 1
Started by Deleted User
Hello,
I have to subset a set of images with the same extentions (coordinates). So, I used the script below which I found on the help. The execution halted at the last instruction with this error message : "Keyword SPATIALREF not allowed in call to: ENVIRASTER::SUBSET".
Could you help me to overcome this error
Thank you.
The script:
Raster = e.OpenRaster(File)
; This is the area of interest:
UpperLeftLat = 35.1
UpperLeftLon = -112.1
LowerRightLat = 34.7
LowerRig...
|
|
|
|
1 |
2314 |
09 Mar 2017 06:02 AM |
|
0 Replies and 2957 Views
the Solution of the equation 2957 0
Started by Deleted User
My newtfunc contain the varion and unknown solution,but the newton method in IDL just allow given initially solution of the equation.
How to get the the Solution of the equation Thanks
Pro TSM_lixiaModel
;http://www.harrisgeospatial.com/docs/ProgrammingGuideTileIterators.html
!Except=0
Compile_opt idl2
e=Envi(/headless)
starttime = Systime(/second)
dir='I:\lwkDATA\waterRegion\MODIS\'
files=file_search(dir,'*RC.tif',count=num)
for i=0,num-1 do begin
file=files&...
|
|
|
|
0 |
2957 |
01 Mar 2017 12:32 AM |
|
5 Replies and 4952 Views
array memory allocation 4952 5
Started by Deleted User
I am reading from a netcdf file with 97 arrays of size (1165,720,1440), reading each array using ncdf_varget. Because of memory allocation problems, I started using delvar to erase the arrays after finished using them - ultimately creating 19 new arrays of the same size from the original arrays (by adding some of them together). Yet am still getting the error Unable to allocate memory: to make array. Cannot allocate memory. When I test delvar and try to print out one of the deleted arrays, th...
|
|
|
|
5 |
4952 |
28 Feb 2017 08:28 AM |
|
2 Replies and 2898 Views
Georeference Meris image with IDL 2898 2
Started by Deleted User
Hi,
I have to write an IDL script to automatise the processing of Meris images. One of the first step is the georeferencement.
In ENVI, it's done by typing this commands :
Geometric correction>Georeference by sensor>Georeference Meris.
But, I don't find the equivalent instructions in IDL.
Does someone have an idea
Thank you.
|
|
|
|
2 |
2898 |
25 Feb 2017 07:51 AM |
|
4 Replies and 2964 Views
rebin 2964 4
Started by Deleted User
I am rebinning from a 0.25 degree lon/lat grid (1440 x 720) to a half degree (720 x 360), using vu_regrid = rebin(vu,720,360,1165, /SAMPLE). In the original, bad (i.e. ocean) values are indicated by 1e20. After rebinning, these values are 2e20. It's supposed to be averaging them - why is it doubling those large values Hard to tell what it's doing with the other values, but just doing quick sum, seems to be ok as the sum is about 4 times less than summing the higher resolution grid.
|
|
|
|
4 |
2964 |
24 Feb 2017 09:08 AM |
|
1 Replies and 2659 Views
warp with tie points 2659 1
Started by Deleted User
Need to warp AWHRR files so that the resulting grid mesh cells borders coincide for all the scenes.
I'm using the following code for warping (in cycle for every file):
; setting projection
units = envi_proj_units_translate('Degrees')
proj = envi_proj_create(/geographic, unit=units)
ENVI_DOIT, 'ENVI_AVHRR_WARP_DOIT', AVHRR_FID=fidin, $
DIMS=dims, FID=fidcal, /IN_MEMORY, BACKGROUND=nodata, $
METHOD=8, OUT_BNAME=&910,1,2,3,4&93, POS=pos, $
PIXEL_SIZE=&910.02, 0.02&93, grid=...
|
|
|
|
1 |
2659 |
23 Feb 2017 04:07 PM |
|
0 Replies and 2524 Views
Only able to store a few days of command history 2524 0
Started by Deleted User
Question:
An IDL user is running into an issue in which he cannot store more than a few days worth of command history in the IDL workbench. He would like to store much more.
Possible answer:
You can increase the amount of history being stored by going to Preferences > IDL, and then changing "History size (Lines)" to a bigger number.
|
|
|
|
0 |
2524 |
23 Feb 2017 01:50 PM |
|
14 Replies and 3862 Views
Removing pixels w/ a value < ... from total histogram pixel count 3862 14
Started by Deleted User
Hello,
In a semi-simplified form here is my question:
Currently a portion of my .pro that outputs the final required value uses the ENVI_STATS_DOIT function to create a histogram of a single band image and then calculates the percent of pixels that fall within the two categories (pixels that meet a condition, pixels that do not) by using hist&911&93/total(hist)*100. Up until not, it worked great and allowed me to determine what percentage of the image was represented by "1's" (oth...
|
|
|
|
14 |
3862 |
17 Feb 2017 01:42 PM |
|
4 Replies and 2662 Views
Contour: difference between IDL 8.0 and IDL 8.4 2662 4
Started by Deleted User
Ciao,
I've upgraded from IDL 8.0 to IDL 8.4, and found a strange difference in the behaviour of CONTOUR.
Here's a simple test:
data = RANDOMU(seed, 9, 9)
smooth = CONTOUR(MIN_CURVE_SURF(data), TITLE='Smoothed', RGB_TABLE=1, /FILL, c_value=&910, 0.2, 0.3, 0.5&93)
If you launch it in IDL 8.0, you get what you (more or less) expect: three filled contours, while the outer part is the white transparent background. I need to include the '0' level, otherwise the inner contour is...
|
|
|
|
4 |
2662 |
16 Feb 2017 06:32 AM |
|
1 Replies and 3259 Views
write a binary file 3259 1
Started by Deleted User
dear all
I have a density file ( array file) with x-dimension=300 and y-demnsion=240,
density=array &91300,240&93.
i want to write it as a binary file with the double precesion format. i used the following commands:
openw, lun, 'density.dat', /get_lun
writeu, lun, density
close, lun
but the problem is that if I open the binary file with the MATLAB the dimension of the density, now, is
density=array &916,12000&93.
I know I have to do some manipulation for writing the binary file to k...
|
|
|
|
1 |
3259 |
14 Feb 2017 09:08 AM |
|