Странные растры в QGIS + GDAL - давайте разберемся!

Сообщение **manisha** » 16 июл 2010, 01:08

Как хорошо известно, пирамидные слои сильно ускоряют отрисовку изображения. Всегда ли? Отнюдь.

Есть например вот такой http://gis-lab.info/data/samples/l7_142-021_1999-187.7z растр. Сцена ландсата. Попробуйте загрузить ее в QGIS (можно убедиться, что пирамиды в ней построены) и поработать с ней. Тормозит?!

В чем причина - в странной структуре растра. Блоки не квадратные как у обычных растров, а размеров ширина растра * 1 пиксель. Смотрим утилитой gdalinfo этот растр и любой другой ландсат, отрисовывающийся быстро. Назовем их условно SLOW и FAST.

SLOW

Код: Выделить всё

Driver: GTiff/GeoTIFF
Files: /data/l7_142-021_1999-187.tif
Size is 8841, 8041
Coordinate System is:
PROJCS["WGS 84 / UTM zone 46N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.2572235629972,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",93],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32646"]]
Origin = (263985.000000000000000,6936015.000000000000000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  263985.000, 6936015.000) ( 88d25'8.31"E, 62d28'47.57"N)
Lower Left  (  263985.000, 6694785.000) ( 88d43'30.72"E, 60d19'14.60"N)
Upper Right (  529215.000, 6936015.000) ( 93d34'5.25"E, 62d33'14.00"N)
Lower Right (  529215.000, 6694785.000) ( 93d31'48.02"E, 60d23'18.23"N)
Center      (  396600.000, 6815400.000) ( 91d 3'35.99"E, 61d27'30.70"N)
Band 1 Block=8841x1 Type=Byte, ColorInterp=Gray
  NoData Value=0
  Overviews: 4421x4021, 2211x2011, 1106x1006, 553x503, 277x252
Band 2 Block=8841x1 Type=Byte, ColorInterp=Undefined
  NoData Value=0
  Overviews: 4421x4021, 2211x2011, 1106x1006, 553x503, 277x252
Band 3 Block=8841x1 Type=Byte, ColorInterp=Undefined
  NoData Value=0
  Overviews: 4421x4021, 2211x2011, 1106x1006, 553x503, 277x252
Band 4 Block=8841x1 Type=Byte, ColorInterp=Undefined
  NoData Value=0
  Overviews: 4421x4021, 2211x2011, 1106x1006, 553x503, 277x252
Band 5 Block=8841x1 Type=Byte, ColorInterp=Undefined
  NoData Value=0
  Overviews: 4421x4021, 2211x2011, 1106x1006, 553x503, 277x252
Band 6 Block=8841x1 Type=Byte, ColorInterp=Undefined
  NoData Value=0
  Overviews: 4421x4021, 2211x2011, 1106x1006, 553x503, 277x252

FAST

Код: Выделить всё

Driver: HFA/Erdas Imagine Images (.img)
Files: /data/LE71590182009205ASN00.img
Size is 8701, 7781
Coordinate System is:
PROJCS["UTM Zone 42, Northern Hemisphere",
    GEOGCS["WGS_1984",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.2572235629972],
            TOWGS84[0,0,0,0,0,0,0]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",69],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["Meter",1],
    AUTHORITY["EPSG","32642"]]
Origin = (423584.999999999941792,6778215.000000000000000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Corner Coordinates:
Upper Left  (  423585.000, 6778215.000) ( 67d34'52.42"E, 61d 7'51.84"N)
Lower Left  (  423585.000, 6544785.000) ( 67d40'6.62"E, 59d 2'8.40"N)
Upper Right (  684615.000, 6778215.000) ( 72d25'28.92"E, 61d 5'42.45"N)
Lower Right (  684615.000, 6544785.000) ( 72d12'51.83"E, 59d 0'9.44"N)
Center      (  554100.000, 6661500.000) ( 69d58'21.00"E, 60d 5'13.26"N)
Band 1 Block=64x64 Type=Byte, ColorInterp=Undefined
  Description = Layer_1
  Min=0.000 Max=159.000 
  Minimum=0.000, Maximum=159.000, Mean=28.300, StdDev=30.435
  Overviews: 2176x1946, 1088x973, 544x487, 272x244, 136x122, 68x61
  Metadata:
    STATISTICS_MINIMUM=0
    STATISTICS_MAXIMUM=159
    STATISTICS_MEAN=28.299991309473
    STATISTICS_MEDIAN=0
    STATISTICS_MODE=0
    STATISTICS_STDDEV=30.435067834072
    STATISTICS_HISTONUMBINS=256
    STATISTICS_HISTOMIN=0
    STATISTICS_HISTOMAX=255
    LAYER_TYPE=athematic
    STATISTICS_HISTOBINVALUES=36150464|0|0|0|128|64|64|0|320|512|0|64|0|0|64|0|192|0|0|64|0|0|0|0|0|0|64|64|0|64|0|0|0|0|0|0|0|64|128|64|0|0|0|0|0|0|0|0|0|64|448|5248|39808|194560|744384|1868800|2862272|2947776|2659520|2677184|2846784|2682880|2312896|2023360|1879872|1663040|1344960|970944|629440|388352|256512|158656|101952|68544|47104|35904|27520|21248|17792|13632|10816|9216|7360|6016|5312|5120|4864|4096|3072|2560|2624|2432|2944|2176|2496|2496|1984|2112|1856|2112|1792|1600|1408|896|1088|1280|704|1216|896|832|1024|1024|1088|896|960|384|832|256|704|384|192|448|832|256|320|256|256|256|192|256|128|384|384|192|192|320|256|256|64|64|128|64|192|192|64|128|192|128|64|0|128|128|128|0|0|0|64|0|0|64|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|
<GDALRasterAttributeTable Row0Min="0" BinSize="1">
  <FieldDefn index="0">
    <Name>Histogram</Name>
    <Type>1</Type>
    <Usage>0</Usage>
  </FieldDefn>
  <Row index="0">
    <F>36150464</F>
  </Row>
  <Row index="1">
    <F>0</F>
  </Row>
  <Row index="2">
    <F>0</F>
  </Row>
  <Row index="3">
    <F>0</F>
  </Row>
  <Row index="4">
    <F>128</F>
  </Row>
  <Row index="5">
    <F>64</F>
  </Row>
  <Row index="6">
    <F>64</F>
  </Row>
  <Row index="7">
    <F>0</F>
  </Row>
  <Row index="8">
    <F>320</F>
  </Row>
  <Row index="9">
    <F>512</F>
  </Row>
  <Row index="10">
    <F>0</F>
  </Row>
  <Row index="11">
    <F>64</F>
  </Row>
  <Row index="12">
    <F>0</F>
  </Row>
  <Row index="13">
    <F>0</F>
  </Row>
  <Row index="14">
    <F>64</F>
  </Row>
  <Row index="15">
    <F>0</F>
  </Row>
  <Row index="16">
    <F>192</F>
  </Row>
  <Row index="17">
    <F>0</F>
  </Row>
  <Row index="18">
    <F>0</F>
  </Row>
  <Row index="19">
    <F>64</F>
  </Row>
  <Row index="20">
    <F>0</F>
  </Row>
  <Row index="21">
    <F>0</F>
  </Row>
  <Row index="22">
    <F>0</F>
  </Row>
  <Row index="23">
    <F>0</F>
  </Row>
  <Row index="24">
    <F>0</F>
  </Row>
  <Row index="25">
    <F>0</F>
  </Row>
  <Row index="26">
    <F>64</F>
  </Row>
  <Row index="27">
    <F>64</F>
  </Row>
  <Row index="28">
    <F>0</F>
  </Row>
  <Row index="29">
    <F>64</F>
  </Row>
  <Row index="30">
    <F>0</F>
  </Row>
  <Row index="31">
    <F>0</F>
  </Row>
  <Row index="32">
    <F>0</F>
  </Row>
  <Row index="33">
    <F>0</F>
  </Row>
  <Row index="34">
    <F>0</F>
  </Row>
  <Row index="35">
    <F>0</F>
  </Row>
  <Row index="36">
    <F>0</F>
  </Row>
  <Row index="37">
    <F>64</F>
  </Row>
  <Row index="38">
    <F>128</F>
  </Row>
  <Row index="39">
    <F>64</F>
  </Row>
  <Row index="40">
    <F>0</F>
  </Row>
  <Row index="41">
    <F>0</F>
  </Row>
  <Row index="42">
    <F>0</F>
  </Row>
  <Row index="43">
    <F>0</F>
  </Row>
  <Row index="44">
    <F>0</F>
  </Row>
  <Row index="45">
    <F>0</F>
  </Row>
  <Row index="46">
    <F>0</F>
  </Row>
  <Row index="47">
    <F>0</F>
  </Row>
  <Row index="48">
    <F>0</F>
  </Row>
  <Row index="49">
    <F>64</F>
  </Row>
  <Row index="50">
    <F>448</F>
  </Row>
  <Row index="51">
    <F>5248</F>
  </Row>
  <Row index="52">
    <F>39808</F>
  </Row>
  <Row index="53">
    <F>194560</F>
  </Row>
  <Row index="54">
    <F>744384</F>
  </Row>
  <Row index="55">
    <F>1868800</F>
  </Row>
  <Row index="56">
    <F>2862272</F>
  </Row>
  <Row index="57">
    <F>2947776</F>
  </Row>
  <Row index="58">
    <F>2659520</F>
  </Row>
  <Row index="59">
    <F>2677184</F>
  </Row>
  <Row index="60">
    <F>2846784</F>
  </Row>
  <Row index="61">
    <F>2682880</F>
  </Row>
  <Row index="62">
    <F>2312896</F>
  </Row>
  <Row index="63">
    <F>2023360</F>
  </Row>
  <Row index="64">
    <F>1879872</F>
  </Row>
  <Row index="65">
    <F>1663040</F>
  </Row>
  <Row index="66">
    <F>1344960</F>
  </Row>
  <Row index="67">
    <F>970944</F>
  </Row>
  <Row index="68">
    <F>629440</F>
  </Row>
  <Row index="69">
    <F>388352</F>
  </Row>
  <Row index="70">
    <F>256512</F>
  </Row>
  <Row index="71">
    <F>158656</F>
  </Row>
  <Row index="72">
    <F>101952</F>
  </Row>
  <Row index="73">
    <F>68544</F>
  </Row>
  <Row index="74">
    <F>47104</F>
  </Row>
  <Row index="75">
    <F>35904</F>
  </Row>
  <Row index="76">
    <F>27520</F>
  </Row>
  <Row index="77">
    <F>21248</F>
  </Row>
  <Row index="78">
    <F>17792</F>
  </Row>
  <Row index="79">
    <F>13632</F>
  </Row>
  <Row index="80">
    <F>10816</F>
  </Row>
  <Row index="81">
    <F>9216</F>
  </Row>
  <Row index="82">
    <F>7360</F>
  </Row>
  <Row index="83">
    <F>6016</F>
  </Row>
  <Row index="84">
    <F>5312</F>
  </Row>
  <Row index="85">
    <F>5120</F>
  </Row>
  <Row index="86">
    <F>4864</F>
  </Row>
  <Row index="87">
    <F>4096</F>
  </Row>
  <Row index="88">
    <F>3072</F>
  </Row>
  <Row index="89">
    <F>2560</F>
  </Row>
  <Row index="90">
    <F>2624</F>
  </Row>
  <Row index="91">
    <F>2432</F>
  </Row>
  <Row index="92">
    <F>2944</F>
  </Row>
  <Row index="93">
    <F>2176</F>
  </Row>
  <Row index="94">
    <F>2496</F>
  </Row>
  <Row index="95">
    <F>2496</F>
  </Row>
  <Row index="96">
    <F>1984</F>
  </Row>
  <Row index="97">
    <F>2112</F>
  </Row>
  <Row index="98">
    <F>1856</F>
  </Row>
  <Row index="99">
    <F>2112</F>
  </Row>
  <Row index="100">
    <F>1792</F>
  </Row>
  <Row index="101">
    <F>1600</F>
  </Row>
  <Row index="102">
    <F>1408</F>
  </Row>
  <Row index="103">
    <F>896</F>
  </Row>
  <Row index="104">
    <F>1088</F>
  </Row>
  <Row index="105">
    <F>1280</F>
  </Row>
  <Row index="106">
    <F>704</F>
  </Row>
  <Row index="107">
    <F>1216</F>
  </Row>
  <Row index="108">
    <F>896</F>
  </Row>
  <Row index="109">
    <F>832</F>
  </Row>
  <Row index="110">
    <F>1024</F>
  </Row>
  <Row index="111">
    <F>1024</F>
  </Row>
  <Row index="112">
    <F>1088</F>
  </Row>
  <Row index="113">
    <F>896</F>
  </Row>
  <Row index="114">
    <F>960</F>
  </Row>
  <Row index="115">
    <F>384</F>
  </Row>
  <Row index="116">
    <F>832</F>
  </Row>
  <Row index="117">
    <F>256</F>
  </Row>
  <Row index="118">
    <F>704</F>
  </Row>
  <Row index="119">
    <F>384</F>
  </Row>
  <Row index="120">
    <F>192</F>
  </Row>
  <Row index="121">
    <F>448</F>
  </Row>
  <Row index="122">
    <F>832</F>
  </Row>
  <Row index="123">
    <F>256</F>
  </Row>
  <Row index="124">
    <F>320</F>
  </Row>
  <Row index="125">
    <F>256</F>
  </Row>
  <Row index="126">
    <F>256</F>
  </Row>
  <Row index="127">
    <F>256</F>
  </Row>
  <Row index="128">
    <F>192</F>
  </Row>
  <Row index="129">
    <F>256</F>
  </Row>
  <Row index="130">
    <F>128</F>
  </Row>
  <Row index="131">
    <F>384</F>
  </Row>
  <Row index="132">
    <F>384</F>
  </Row>
  <Row index="133">
    <F>192</F>
  </Row>
  <Row index="134">
    <F>192</F>
  </Row>
  <Row index="135">
    <F>320</F>
  </Row>
  <Row index="136">
    <F>256</F>
  </Row>
  <Row index="137">
    <F>256</F>
  </Row>
  <Row index="138">
    <F>64</F>
  </Row>
  <Row index="139">
    <F>64</F>
  </Row>
  <Row index="140">
    <F>128</F>
  </Row>
  <Row index="141">
    <F>64</F>
  </Row>
  <Row index="142">
    <F>192</F>
  </Row>
  <Row index="143">
    <F>192</F>
  </Row>
  <Row index="144">
    <F>64</F>
  </Row>
  <Row index="145">
    <F>128</F>
  </Row>
  <Row index="146">
    <F>192</F>
  </Row>
  <Row index="147">
    <F>128</F>
  </Row>
  <Row index="148">
    <F>64</F>
  </Row>
  <Row index="149">
    <F>0</F>
  </Row>
  <Row index="150">
    <F>128</F>
  </Row>
  <Row index="151">
    <F>128</F>
  </Row>
  <Row index="152">
    <F>128</F>
  </Row>
  <Row index="153">
    <F>0</F>
  </Row>
  <Row index="154">
    <F>0</F>
  </Row>
  <Row index="155">
    <F>0</F>
  </Row>
  <Row index="156">
    <F>64</F>
  </Row>
  <Row index="157">
    <F>0</F>
  </Row>
  <Row index="158">
    <F>0</F>
  </Row>
  <Row index="159">
    <F>64</F>
  </Row>
  <Row index="160">
    <F>0</F>
  </Row>
  <Row index="161">
    <F>0</F>
  </Row>
  <Row index="162">
    <F>0</F>
  </Row>
  <Row index="163">
    <F>0</F>
  </Row>
  <Row index="164">
    <F>0</F>
  </Row>
  <Row index="165">
    <F>0</F>
  </Row>
  <Row index="166">
    <F>0</F>
  </Row>
  <Row index="167">
    <F>0</F>
  </Row>
  <Row index="168">
    <F>0</F>
  </Row>
  <Row index="169">
    <F>0</F>
  </Row>
  <Row index="170">
    <F>0</F>
  </Row>
  <Row index="171">
    <F>0</F>
  </Row>
  <Row index="172">
    <F>0</F>
  </Row>
  <Row index="173">
    <F>0</F>
  </Row>
  <Row index="174">
    <F>0</F>
  </Row>
  <Row index="175">
    <F>0</F>
  </Row>
  <Row index="176">
    <F>0</F>
  </Row>
  <Row index="177">
    <F>0</F>
  </Row>
  <Row index="178">
    <F>0</F>
  </Row>
  <Row index="179">
    <F>0</F>
  </Row>
  <Row index="180">
    <F>0</F>
  </Row>
  <Row index="181">
    <F>0</F>
  </Row>
  <Row index="182">
    <F>0</F>
  </Row>
  <Row index="183">
    <F>0</F>
  </Row>
  <Row index="184">
    <F>0</F>
  </Row>
  <Row index="185">
    <F>0</F>
  </Row>
  <Row index="186">
    <F>0</F>
  </Row>
  <Row index="187">
    <F>0</F>
  </Row>
  <Row index="188">
    <F>0</F>
  </Row>
  <Row index="189">
    <F>0</F>
  </Row>
  <Row index="190">
    <F>0</F>
  </Row>
  <Row index="191">
    <F>0</F>
  </Row>
  <Row index="192">
    <F>0</F>
  </Row>
  <Row index="193">
    <F>0</F>
  </Row>
  <Row index="194">
    <F>0</F>
  </Row>
  <Row index="195">
    <F>0</F>
  </Row>
  <Row index="196">
    <F>0</F>
  </Row>
  <Row index="197">
    <F>0</F>
  </Row>
  <Row index="198">
    <F>0</F>
  </Row>
  <Row index="199">
    <F>0</F>
  </Row>
  <Row index="200">
    <F>0</F>
  </Row>
  <Row index="201">
    <F>0</F>
  </Row>
  <Row index="202">
    <F>0</F>
  </Row>
  <Row index="203">
    <F>0</F>
  </Row>
  <Row index="204">
    <F>0</F>
  </Row>
  <Row index="205">
    <F>0</F>
  </Row>
  <Row index="206">
    <F>0</F>
  </Row>
  <Row index="207">
    <F>0</F>
  </Row>
  <Row index="208">
    <F>0</F>
  </Row>
  <Row index="209">
    <F>0</F>
  </Row>
  <Row index="210">
    <F>0</F>
  </Row>
  <Row index="211">
    <F>0</F>
  </Row>
  <Row index="212">
    <F>0</F>
  </Row>
  <Row index="213">
    <F>0</F>
  </Row>
  <Row index="214">
    <F>0</F>
  </Row>
  <Row index="215">
    <F>0</F>
  </Row>
  <Row index="216">
    <F>0</F>
  </Row>
  <Row index="217">
    <F>0</F>
  </Row>
  <Row index="218">
    <F>0</F>
  </Row>
  <Row index="219">
    <F>0</F>
  </Row>
  <Row index="220">
    <F>0</F>
  </Row>
  <Row index="221">
    <F>0</F>
  </Row>
  <Row index="222">
    <F>0</F>
  </Row>
  <Row index="223">
    <F>0</F>
  </Row>
  <Row index="224">
    <F>0</F>
  </Row>
  <Row index="225">
    <F>0</F>
  </Row>
  <Row index="226">
    <F>0</F>
  </Row>
  <Row index="227">
    <F>0</F>
  </Row>
  <Row index="228">
    <F>0</F>
  </Row>
  <Row index="229">
    <F>0</F>
  </Row>
  <Row index="230">
    <F>0</F>
  </Row>
  <Row index="231">
    <F>0</F>
  </Row>
  <Row index="232">
    <F>0</F>
  </Row>
  <Row index="233">
    <F>0</F>
  </Row>
  <Row index="234">
    <F>0</F>
  </Row>
  <Row index="235">
    <F>0</F>
  </Row>
  <Row index="236">
    <F>0</F>
  </Row>
  <Row index="237">
    <F>0</F>
  </Row>
  <Row index="238">
    <F>0</F>
  </Row>
  <Row index="239">
    <F>0</F>
  </Row>
  <Row index="240">
    <F>0</F>
  </Row>
  <Row index="241">
    <F>0</F>
  </Row>
  <Row index="242">
    <F>0</F>
  </Row>
  <Row index="243">
    <F>0</F>
  </Row>
  <Row index="244">
    <F>0</F>
  </Row>
  <Row index="245">
    <F>0</F>
  </Row>
  <Row index="246">
    <F>0</F>
  </Row>
  <Row index="247">
    <F>0</F>
  </Row>
  <Row index="248">
    <F>0</F>
  </Row>
  <Row index="249">
    <F>0</F>
  </Row>
  <Row index="250">
    <F>0</F>
  </Row>
  <Row index="251">
    <F>0</F>
  </Row>
  <Row index="252">
    <F>0</F>
  </Row>
  <Row index="253">
    <F>0</F>
  </Row>
  <Row index="254">
    <F>0</F>
  </Row>
  <Row index="255">
    <F>0</F>
  </Row>
</GDALRasterAttributeTable>

и далее аналогично по оставшимся каналам...

Видно, что в первом случае блоки Block=8841x1, а во втором Block=64x64.

Теперь собственно вопросы:

Как и почему они такие получились? Исходники - поканальные тифы из гловиса, склеивал или через gdaltools, или эрдасом, не помню точно. Впрочем и тем. и другим склеивал не раз - не было такого.

Как переконвертировать к нормальной структуре? Пересохранение в другой формат (tif -> img и наоборот) не помогает, перегонка в другую проекцию тоже.

Можно ли ускорить отображение таких растров?

Интересно также быстро ли их отрисовывают другие программы. Помниться, и раньше попадались некоторые растры, которые очень медленно отрисовывались в трешке и девятке, не вэтом ли причина?

Сообщение **manisha** » 16 июл 2010, 01:30

Итак, выяснили сразу две вещи:

1) растры со странной структурой беруться в результате смены проекции: смены 46-ой зоны UTM на 46-ю зону UTM. Входит и выходит, как говориться...

2) растр не тормозит :) тормозят лю.. тьфу, векторы. В этом же проекте есть векторный слой, с границами субъектов. Я как то даже забыл про него, так как он смехотворно мал и в других программах никак не сказывается на скорости отображения. В символики слой без заливки, только с границей. Так вот, на крупных масштабах - начиная примерно с 1:200 000 и крупнее - начинаются феноменальные тормоза. Причина видимо в том, что слой на всю РФ, а проекция - UTM...

Максим Дубинин

manisha писал(а):смены 46-ой зоны UTM на 46-ю зону UTM

а что это за смена такая?

Форумы GIS-Lab.info

Странные растры в QGIS + GDAL - давайте разберемся!

Странные растры в QGIS + GDAL - давайте разберемся!

Re: Странные растры в QGIS + GDAL - давайте разберемся!

Re: Странные растры в QGIS + GDAL - давайте разберемся!

Кто сейчас на конференции