I want to make a script that calculates certain characteristics in a satellite image for the day April 3, 2004.
The first thing I want to do:
1. how to select the Region of Mexico from a script?
How to convert lines and elements of an image to coordinates of latitude and longitude?
There is some method or function to select the domain mentioned?
2. There is some method or function to display this data from script (for example: hourly data for the day April 3 2004) on screen of McIDAS-V?
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?
4. What is the best way to add the features mentioned at the code below
Code: Select all
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();
# now look through the first band y count pixels
# MCS detected when his temperature infrared (TIR) is < 219 K
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")