Scribd
Upload
45 views

Auto/Semi-Automated Processing of ENVISAT ASAR Imagery of The Arctic Ocean

This document discusses automating the processing of ENVISAT ASAR Arctic Ocean imagery into ArcGIS. It involves: 1. Downloading imagery from the PolarView archive and uncompressing the files. 2. Converting the ASAR images to rasters, extracting the image boundaries as polygons, and clipping the rasters. 3. Defining the correct projection and placing the rasters and polylines into a file geodatabase. The process is demonstrated on sample ENVISAT imagery from a 2010 Chinese Arctic expedition, and estimates are provided for processing typical daily Arctic Ocean datasets.

Uploaded by

Buuvee Be
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
45 views

Auto/Semi-Automated Processing of ENVISAT ASAR Imagery of The Arctic Ocean

This document discusses automating the processing of ENVISAT ASAR Arctic Ocean imagery into ArcGIS. It involves: 1. Downloading imagery from the PolarView archive and uncompressing the files. 2. Converting the ASAR images to rasters, extracting the image boundaries as polygons, and clipping the rasters. 3. Defining the correct projection and placing the rasters and polylines into a file geodatabase. The process is demonstrated on sample ENVISAT imagery from a 2010 Chinese Arctic expedition, and estimates are provided for processing typical daily Arctic Ocean datasets.

Uploaded by

Buuvee Be
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 41

Auto/Semi-Automated processing of ENVISAT ASAR imagery of the Arctic Ocean

- Explore ArcGIS manual procedure

- Implement as a Python script

XUELONG cruise during CHINARE-2010 from July 1 Xianmen to Sept 20 Shanghai

http://www.utsa.edu/LRSG/Arctic/index.htm

XUELONG, crew, and staff at Long-term Ice Station, Aug 17, 2010

http://www.utsa.edu/LRSG/Arctic/index.htm

"Summer sea ice motion..." by Ron Kwok, JPL Geophysical Research Letters, Vol 35, 2008
Jul Aug

2003

Note: Longitude lines added for 0, 90, 150, 180, and 270 degrees Purple lines mark the main region covered during by the XUELONG during CHINARE-2010

2004

"Summer sea ice motion..." by Ron Kwok, JPL Geophysical Research Letters, Vol 35, 2008
Jul Aug

2005

Note: Longitude lines added for 0, 90, 150, 180, and 270 degrees Purple lines mark the main region covered during by the XUELONG during CHINARE-2010

2006

Background for CHINARE-2010


While the XUELONG icebreaker went by sea, ENVISAT took SAR imagery from space XUELONG saw sea ice from Jul 21 to Aug 28 160 ENVISAT ASAR images covering Jul 23 to Sep 28, plus 6 missing or damaged images Data for one image takes up 100 to 500 MB

ENVISAT drawing

http://darc.nerc.ac.uk/asset/images/envisat.jpg

Desired results
Automated process to create image rasters and boundary polylines and put them into a geodatabase Reasonable to handle ENVISAT ASAR datasets as fast as they are produced by the satellite Average for Beaufort Sea (Arctic Ocean) only 16 daily datasets ~ 0.8 GB compressed data

Average for entire Arctic and Antarctic Oceans 80 daily datasets ~ 4 GB compressed data

Step 1
Use web browser Get ENVISAT ASAR data

In class: one dataset


Homework: multiple datasets

>>> www.polarview.aq <<< Current ENVISAT ASAR datasets

PolarView holds data only 1 month Can filter by region

Click an image to preview

Can display list of images Avg 16 daily datasets for Beaufort Sea

Arctic Ocean - sample 24 hours 35 datasets

Antarctic Ocean - sample 24 hours 28 datasets

Select and download datasets

Step 2
Uncompress and display original files

ESA compressed the datafiles twice

Uncompress twice

Open the nested folders

Original ASAR images are very dark

Fix brightness/contrast

All original files opened in ArcGIS

Step 3
Convert ASAR image area to raster of a polygon

Raster image pixel value: 16 bit range Background/radar shadow data value = 0

Reclassify Tool (Spatial Analyst)

Raster polygon data value = 0 Background/shadow data value = NoData

Step 4
Convert raster polygon to a shapefile polygon

Raster to Polygon Tool (Conversion) (Note: Simply polygons off)

Shapefile polygon data value = 0 Background/shadow data value = NoData

Step 5
Convert shapefile polygon to shapefile polyline

Feature to Line (Data Management)

Shapefile polyline created

Files after creating image outline

Step 6
Clip background, keep only ASAR image data

Clip Tool (Data Management)

Background/shadow data value = NoData All other pixels = original ASAR image

ArcMap - clipped image and polyline (Note: No spatial reference)

Define Projection Tool (Data Management)

What happens if you select a different polar projection?

ArcMap - clipped image and polyline (With correct spatial reference)

Sample running of the script


Import required modules Tar-zip inbox folder Tar-zip done folder Working data folder Geodatabase outputs Coordinate system F:/Chinare_2010_ASAR/TarZipInbox F:/Chinare_2010_ASAR/TarZipDone F:/Chinare_2010_ASAR/Working F:/Chinare_2010_ASAR/Arctic_ArcGIS_10.gdb WGS 1984 NSIDC Sea Ice Polar Stereographic North.prj

Start processing loop


Working on dataset 1 of 1 named WSM_SS_20110420_220800_4276_4 Extract from tarzip file to dataset folder Move ASAR data up to the dataset folder Original spatial reference - Unknown Redefine as - WGS 1984 NSIDC Sea Ice Polar Stereographic North.prj Calculate statistics Build pyramids Reclassify 0 to NoData and positive to 0 Convert reclassified raster to polygon Convert polygon to polyline Use polygon to clip original raster Build pyramids for clipped raster Copy polyline to output geodatabase F:/Chinare_2010_ASAR/Working/WSM_SS_20110420_220800_4276_4\WSM_SS_20110420_22080 0_4276_4_polyline.shp Successfully converted: F:/Chinare_2010_ASAR/Arctic_ArcGI S_10.gdb\WSM_SS_20110420_220800_4276_4_polyline Copy clipped raster to output geodatabase Successfully converted: F:/Chinare_2010_ASAR/Working/WSM_SS_20110420_220800_427 6_4\WSM_SS_20110420_220800_4276_4_clip.tif To F:/Chinare_2010_ASAR/Arctic_ArcGIS _10.gdb\WSM_SS_20110420_220800_4276_4_clip Finish processing loop

Estimates for Beaufort Sea (Arctic Ocean) Average 16 daily datasets on PolarView
Downloading time not included
Processing time 80 minutes Mac laptop with 1.8 GHz processor, running Windows 7 All data stored on external SSD via USB 2.0 connection Approximate data sizes Compressed dataset from PolarView Uncompressed ASAR datafile (GeoTIFF) All ArcGIS working files after processing Image boundary and clipped raster MB 800 2320 3360 960

Questions?

You might also like