欧美bbbwbbbw肥妇,免费乱码人妻系列日韩,一级黄片

Python中ArcPy柵格裁剪柵格(批量對齊柵格圖像范圍并統(tǒng)一行數(shù)與列數(shù))

 更新時(shí)間:2023年02月24日 14:26:03   作者:瘋狂學(xué)習(xí)GIS  
本文介紹基于Python中ArcPy模塊,實(shí)現(xiàn)基于柵格圖像批量裁剪柵格圖像,同時(shí)對齊各個(gè)柵格圖像的空間范圍,統(tǒng)一其各自行數(shù)與列數(shù)的方法,本文通過實(shí)例代碼給大家介紹的非常詳細(xì),需要的朋友可以參考下

  本文介紹基于PythonArcPy模塊,實(shí)現(xiàn)基于柵格圖像批量裁剪柵格圖像,同時(shí)對齊各個(gè)柵格圖像的空間范圍統(tǒng)一其各自行數(shù)列數(shù)的方法。

  首先明確一下我們的需求?,F(xiàn)有某一地區(qū)的多張柵格遙感影像,其雖然都大致對應(yīng)著同樣的地物范圍,但不同柵格影像之間的空間范圍、行數(shù)列數(shù)、像元的位置等都不完全一致;例如,某一景柵格影像會比其他柵格影像多出一行,而另一景柵格影像可能又會比其他柵格影像少一列等等。我們希望可以以其中某一景柵格影像為標(biāo)準(zhǔn),將全部的柵格影像的具體范圍、行數(shù)、列數(shù)等加以統(tǒng)一。

  本文所用到的具體代碼如下。

# -*- coding: utf-8 -*-
"""
Created on Thu Dec 29 21:13:19 2022

@author: fkxxgis
"""

import arcpy

tif_file_path = r"E:\02_Project\01_Chlorophyll\ClimateZone\Original"
result_file_path = r"E:\02_Project\01_Chlorophyll\ClimateZone\Original_Snap/"
snap_file_name = r"E:\02_Project\01_Chlorophyll\ClimateZone\Original\F_LC.tif"

arcpy.env.workspace = tif_file_path
arcpy.env.snapRaster = snap_file_name

tif_file_list = arcpy.ListRasters("*", "tif")

for tif_file in tif_file_list:
    key_name = tif_file.split(".tif")[0] + "S.tif"
    arcpy.Clip_management(tif_file,
                          "#",
                          result_file_path + key_name,
                          snap_file_name,
                          "#",
                          "#",
                          "MAINTAIN_EXTENT")

  其中,tif_file_path是保存有我們原有柵格圖像的路徑,result_file_path是裁剪后各個(gè)結(jié)果圖像的保存路徑(記得在這一路徑后加一個(gè)正斜杠/,否則之后輸出結(jié)果的路徑會有問題),snap_file_name是裁剪其他柵格圖像時(shí),所用的模板柵格圖像——因?yàn)槲覀円y(tǒng)一各個(gè)柵格圖像的行號與列號,所以很顯然,這里這個(gè)模板圖像就需要找各個(gè)柵格圖像中,行數(shù)與列數(shù)均為最少的那一景圖像。這里需要注意,如果大家的各個(gè)柵格圖像中,行數(shù)與列數(shù)最少的柵格不是同一個(gè)柵格,那么可以分別用行數(shù)最少、列數(shù)最少的這兩個(gè)柵格分別作為模板,執(zhí)行兩次上述代碼。

  代碼整體思路也很簡單:首先,我們基于arcpy.ListRasters()函數(shù),獲取tif_file_path路徑下原有的全部.tif格式的圖像文件,并以列表的形式存放于tif_file_list中;隨后,逐一取出tif_file_list列表中的柵格文件,進(jìn)行裁剪處理。這里的裁剪我們是通過arcpy.Clip_management()函數(shù)來實(shí)現(xiàn)的,其各項(xiàng)參數(shù)的具體含義大家可以參考官方幫助文檔,我們這里就只對本文中需要修改的參數(shù)加以介紹。

  其中,第一個(gè)參數(shù)就是當(dāng)前循環(huán)所用的柵格圖像文件,第三個(gè)參數(shù)是結(jié)果文件的保存路徑與文件名,第四個(gè)參數(shù)則是模板文件;最后一個(gè)參數(shù)"MAINTAIN_EXTENT"是為了保證得到的裁剪后結(jié)果圖像嚴(yán)格與模板圖像的行數(shù)、列數(shù)相匹配。除此之外,幾個(gè)"#"表示我們對其他參數(shù)暫時(shí)不配置。

  此外,在代碼開頭的這句arcpy.env.snapRaster = snap_file_name,表明我們將以所選用的模板文件為標(biāo)準(zhǔn),使得輸出的結(jié)果文件的像元大小、圖像范圍等與模板文件保持一致。這里需要注意,這一句代碼與前述的"MAINTAIN_EXTENT"參數(shù)缺一不可——只有二者同時(shí)出現(xiàn),才可以保證輸出結(jié)果與模板文件是嚴(yán)格一致的。

  另一方面,由于我們用到了ArcPy模塊,因此如果大家的Python版本是3.0及以上,則需要在ArcMap軟件中的Python運(yùn)行框,或其對應(yīng)的IDLE(如下圖所示)中運(yùn)行上述代碼。

  運(yùn)行結(jié)果后,可以發(fā)現(xiàn)所有輸出結(jié)果文件就具有完全一致的行數(shù)與列數(shù)了,且其各自的像元位置也是完全一致的。

到此這篇關(guān)于Python中ArcPy柵格裁剪柵格(批量對齊柵格圖像范圍并統(tǒng)一行數(shù)與列數(shù))的文章就介紹到這了,更多相關(guān)Python arcpy柵格內(nèi)容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家!

相關(guān)文章

  • Pytest如何使用skip跳過執(zhí)行測試

    Pytest如何使用skip跳過執(zhí)行測試

    這篇文章主要介紹了Pytest如何使用skip跳過執(zhí)行測試,文中通過示例代碼介紹的非常詳細(xì),對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下
    2020-08-08
  • 在Pycharm中對代碼進(jìn)行注釋和縮進(jìn)的方法詳解

    在Pycharm中對代碼進(jìn)行注釋和縮進(jìn)的方法詳解

    今天小編就為大家分享一篇在Pycharm中對代碼進(jìn)行注釋和縮進(jìn)的方法詳解,具有很好的參考價(jià)值,希望對大家有所幫助。一起跟隨小編過來看看吧
    2019-01-01
  • Python 中的 import 機(jī)制之實(shí)現(xiàn)遠(yuǎn)程導(dǎo)入模塊

    Python 中的 import 機(jī)制之實(shí)現(xiàn)遠(yuǎn)程導(dǎo)入模塊

    模塊導(dǎo)入( import ),是指在一個(gè)模塊中使用另一個(gè)模塊的代碼的操作,它有利于代碼的復(fù)用。這篇文章主要介紹了Python 的 import 機(jī)制:實(shí)現(xiàn)遠(yuǎn)程導(dǎo)入模塊,需要的朋友可以參考下
    2019-10-10
  • Python OpenCV 使用滑動(dòng)條來調(diào)整函數(shù)參數(shù)的方法

    Python OpenCV 使用滑動(dòng)條來調(diào)整函數(shù)參數(shù)的方法

    這篇文章主要介紹了Python OpenCV 使用滑動(dòng)條來調(diào)整函數(shù)參數(shù)的方法,文中通過示例代碼介紹的非常詳細(xì),對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧
    2019-07-07
  • Pandas缺失值2種處理方式代碼實(shí)例

    Pandas缺失值2種處理方式代碼實(shí)例

    這篇文章主要介紹了Pandas缺失值2種處理方式代碼實(shí)例,文中通過示例代碼介紹的非常詳細(xì),對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下
    2020-06-06
  • python Gabor濾波器講解

    python Gabor濾波器講解

    這篇文章主要介紹了python Gabor濾波器講解,本文給大家介紹的非常詳細(xì),對大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友可以參考下
    2020-10-10
  • 利用Python實(shí)現(xiàn)網(wǎng)絡(luò)測試的示例代碼

    利用Python實(shí)現(xiàn)網(wǎng)絡(luò)測試的示例代碼

    Speedtest CLI 為命令行帶來 Speedtest 背后的可信技術(shù)和全球服務(wù)器網(wǎng)絡(luò)。本文將利用它進(jìn)行網(wǎng)絡(luò)測試,感興趣的小伙伴可以了解一下
    2022-04-04
  • Python中執(zhí)行分位數(shù)回歸的示例詳解

    Python中執(zhí)行分位數(shù)回歸的示例詳解

    分位數(shù)回歸是線性回歸的擴(kuò)展版本,分位數(shù)回歸構(gòu)建一組變量(也稱為自變量)和分位數(shù)(也稱為因變量)之間的關(guān)系,下面我們就來看看Python如何執(zhí)行分位數(shù)回歸吧
    2024-03-03
  • python面試題Python2.x和Python3.x的區(qū)別

    python面試題Python2.x和Python3.x的區(qū)別

    這篇文章主要介紹了python面試題Python2.x和Python3.x的區(qū)別 ,在面試中也經(jīng)常會問到,本文給大家介紹的非常詳細(xì),具有一定的參考借鑒價(jià)值,需要的朋友可以參考下
    2019-05-05
  • python super函數(shù)使用方法詳解

    python super函數(shù)使用方法詳解

    這篇文章主要介紹了python super函數(shù)使用方法詳解,需要的朋友可以參考下
    2020-02-02

最新評論