1 / 20

Sea Ice Maps from NPI

Sea Ice Maps from NPI. Maps produced. 30-year Maximum (for every month ) Areas where at least one day ice last 30 years 30 year Minimum (for every month ) Areas where Every single day ice last 30 years Ice Persistency ( Isfrekvens )

Télécharger la présentation

Sea Ice Maps from NPI

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Sea Ice Maps from NPI

  2. Mapsproduced • 30-year Maximum (for everymonth) • Areas where at leastonedayice last 30 years • 30 year Minimum (for everymonth) • Areas whereEvery single dayice last 30 years • Ice Persistency (Isfrekvens) • Percentageof Days havingice last 30 years

  3. Output in TwoProjections • NSIDC EPSG:3411 and North Pole LAEA Europe EPSG:3575

  4. Basedon NSIDC data • Sea Ice Concentrations from Nimbus-7 SMMR and DMSP SSM/I-SSMIS Passive Microwave Data • 25 km • 26 October 1978 – most current processing • Daily

  5. Filtering NSIDC Data • Original NSIDC data containserroneouspixels at coast line • At coast, usuallyvalues less than 15% ice • Still, thesepixelscauseerror in long-term maximumicemap • Solution: Filteringbasedon area

  6. Filtering NSIDC Data • The following is applied to everydaily NSIDC file: • Step 1: Whole area: Pixelis consideredicewhen present day is ice AND daybefore or after is ice. • Step 2: In addition, Norwegian coast (defined by mask) onlyiceif present dayice and twodaysbefore and afterice as well • Step 3: The Baltic and somecoastal areas iceifpresent dayice and twodaysbefore and afterice as well • Step 4: Areas between UK and Lofoten, inland lakes and Pacific never considered to be ice • This mask can be adjusted to filter some areas in a specialway

  7. FilteringResult – Consecutive Days • Purple area filteredawaywhenonlyiceconsideredifprevious or nextdayalsoice

  8. FilteringResult • Results shows howthe pink maximum line excludescoastal areas and noise

  9. FilteringResult • No Filtering at all versus filtering … biggestdifference in October • Many non-consecutiveicepixels in boundaryzone

  10. Create Max / Min Maps • Create a NumberOfDay Raster • Loop through all files and add «1» whenvalue > 38 ( 15.2%) • Maximum Map shows where at leastonedayice in wholeperiod: • No icewhereNumberOfDay - Raster is zero, otherwiseice • Minimum Map shows whereevery single day in periodice: • Ice whereNumberOfDay – Raster valueequals total numberofdays-

  11. Create Ice Persistency • In everypixelthepercentageofdayswithice cover > 15% • Calculation: • Loop througheach file • If ice in pixel, add 100/NumberOfDays to persistence raster

  12. Filter Noise in Max/Min maps • Remove singular pixelsusinggdal_sieve filter • Possibility to adjust parameters! • Patches in front ofmaximumextent is removed by choosingthelargest polygon • Baltic is alsoremoved, butadded and treatedseparately

  13. Convert Polygon to Line • Convert polygon to line and subtractcoastline to have iceedgeonly

  14. Smooth Lines • QGIS GeneralizerPlug-In with snakes filter givesgood line-fit • Beziersmoothing in ArcGISremovesanyedges

  15. Smoothed Lines • Comparison to published and previouslycalculated (Gunnar’sMatlab script)

  16. Ice PersistencyCompared • ComparingGunnar’sMatlab and Python script

  17. Baltic Area

  18. Alaska Japan China • East Chinese Sea / Yellow Sea has ice in persistence and maximummap from ouralgorithm • NSIDC mapindicatesthatthis is «real»

  19. Python Script • Availableat https://github.com/npolar/RemoteSensing/tree/master/IcePersistency • Script is documentedwithinthecode • Steps: • 1) Run Script • 2) Manually: Apply snakes filtering QGIS • 3) Manually: ApplyBeziersmoothingArcGIS

More Related