This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.
Hi Boris, re: > I want to make a script that calculates certain characteristics in a > satellite image for the day April 3, 2004. The data is downloaded from > Unidata in the following format: 200404031215.goesi12r and are hourly. OK, those images are in McIDAS AREA format. re: > The first thing I want to do: > > 1. how to select the Region of Mexico from a script? Selecting a region of an image for manipulation/analysis is very easy in McIDAS-X. I am not sure of how to do it in McIDAS-V/IDV. Questions: - what calculation(s) are you wanting to do? I see the code snippit you included in your email, but I am not sure if I understand what it is you are trying to do. - do you just want to list out the values for pixels that match certain conditions, or do you want to highlight those pixels in the image display? If you just want to highlight the pixels that meet certain conditions (e.g., values), you could do this by manipulating the color enhancement you are using. re: > Doing some deployments in McIDAS-V and after having researched the types > of coordinates in McIDAS-V , I note that the values ??of lines and elements > which correspond to Region of Mexico are the values ??determined by the > "AREA Coordinates" (see 3abril.jpg in attached), they are the values > ??associated > with the image displayed on screen ? McIDAS (McIDAS-X and McIDAS-V and the IDV) use the coordinate concepts introduced in McIDAS: image coordinates - lines and elements from the original satellite scan AREA coordinates - that portion of a satellite scan that is contained in an AREA file. An AREA file can contain all or part of a satellite scan AREA coordinates - Earth reference for each pixel. When the pixel is out in space, there is no Earth reference Any pixel in an image can be referred to uniquely/correctly in at least two of these coordinate systems. Why 2 and not 3? A portion of a satellite image's scan may be out in space, not over ground. re: > if that's true, How to convert lines and elements of an image to coordinates > of latitude and longitude? This is done by knowing the navigation (projection) transforms. This can be quite complicated for satellite images in their native projection. The IDV/McIDAS-V does this for you on the fly: the equations that represent the navigation transforms are part of the packages. Hmm... if you remapped the images you got into a Lat,Lon projection (easiest done in McIDAS-X). If this is done, the mapping from line,element to lat,lon is simple. re: > There is some method or function to select the domain mentioned? I don't know. I will have to ask our IDV developers. What did the folks at SSEC say? re: > 2. There is some method or function to display this data on screen of > McIDAS-V? I can't speak about McIDAS-V since it is under the control of SSEC, but I would think that the answer is yes. re: > 3. I want to display a image with the values ??of brightness, where the range > is from 201 to 207 and segment the region corresponding to these values. > Is there any way to do this? I don't know; I will have to ask the IDV guys. re: > The script is as follows: > > def scm(directory): > import os; > fs = os.listdir(directory); > > # now fs will be a list of all the files in directory > > from edu.wisc.ssec.mcidas import AreaFile; > for name in fs: > > print "Reading in:",name > af = AreaFile(directory+"/"+name); > ad = af.getAreaDirectory(); > count = 0; > data = af.getFloatData(); > > > # Mesoscale Convective System detected > > for i in xrange(ad.getLines()): > for j in xrange(ad.getElements()): > if 552 < i < 900 and 244 < j < 572: > if (data[0][i][j]) > 206 and (data[0][i][j]) < 208: > print i, j, data[0][i][j]; > > > scm("/home/mcidas/Documents/datos/IR") I understand what this script is supposed to do, but I don't know the basis for it. Can you elaborate? Cheers, Tom -- **************************************************************************** Unidata User Support UCAR Unidata Program (303) 497-8642 P.O. Box 3000 address@hidden Boulder, CO 80307 ---------------------------------------------------------------------------- Unidata HomePage http://www.unidata.ucar.edu **************************************************************************** Ticket Details =================== Ticket ID: NSZ-383912 Department: Support McIDAS Priority: Normal Status: Closed