c++ - Python スクリプト内で gdal_calc を複数回使用する、Visual C++ for Python 2.7 の問題

Question

以下に投稿したスクリプトに関して、現在できることを使い果たしました。このスクリプトで gdal_calc.py を正しく呼び出すことができません。win32 エラーに対する .dll ファイルの不足に関するさまざまなエラー メッセージを受け取りました。

GDAL をインストールして再インストールしましたが、osgeo がないことを示すエラー メッセージが引き続き表示されます。また、Visual C++ for Python に問題があるため、Python27 パスウェイ内に GDAL を適切にインストールできません。適切な .bat ファイルまたはこのメッセージが表示されます。

致命的なエラー C1083: インクルード ファイルを開けません: 'cpl_port.h': そのようなファイルまたはディレクトリはありません

私は初心者であり、これは私の頭の中にあるので、指し示していただけると幸いです。

ありがとう、

import os
import sys
from osgeo import gdal
from gdalconst import *
import numpy as np
import math
import subprocess
from subprocess import call
import math
import PIL
from PIL import Image

# constants

gdalTranslate = r'C:\\Program Files (x86)\\GDAL\\gdal_translate.exe'
src = r"C:\Users\jack\Desktop\RUA_FIL\IMG-HV-ALOS2110871010-160611-HBQR1.5RUA.img"
srcVH = r"C:\Users\jack\Desktop\RUA_FIL\IMG-VH-ALOS2110871010-160611-HBQR1.5RUA.img"
dstVH = r"C:\Users\jack\Desktop\New_Trash\vhFloat32-3.img"
dst = r"C:\\Users\jack\Desktop\\New_Trash\hvFloat32-3.img"
cmd = "-ot Float32 -of HFA" # hopefully this works
#gdalCalc = r'C:\Python27\ArcGIS10.4\Scripts\dist\gdal_calc.exe'
#hvFloat32 = r"C:\Users\jack\Desktop\New_Trash\hvFloat32-1.img"
#vhFloat32 = r"C:\Users\jack\Desktop\New_Trash\hvFloat32-1.img"


#prevName = r"C:\Users\jack\Desktop\RUA_FIL\IMG-HV-ALOS2110871010-160611-HBQR1.5RUA"
#newName = r"C:\Users\jack\Desktop\RUA_FIL\IMG-HV-ALOS2110871010-160611-HBQR1.5RUA.img"

#os.rename(prevName,newName)


# setting def

#def gdal_Translate (item):
    #return item
hvfullCmd = ' '.join([gdalTranslate, cmd, gdal_Translate(src), gdal_Translate(dst)])
subprocess.Popen(hvfullCmd)

vhfullCmd = ' '.join([gdalTranslate,cmd,gdal_Translate(srcVH),gdal_Translate(dstVH)])
subprocess.Popen(vhfullCmd)


# Run gdal_calc.py for HV

try:
    os.system('C:\Users\jack\Python\QGIS Python Scripts\ridofZs.py')
except:
    print ("This isn't working")
    sys.exit()


subprocess.call([sys.executable,'C:\Program Files (x86)\GDAL\gdal_calc.py', '-A', 'C:\Users\jack\Desktop\New_Trash\hvFloat32.img', '--outfile=C:\Users\jack\Desktop\New_Trash\NEWCALCHV-2.img','--calc=A*(A>=1)'])



# Run gdal_calc.py for VH

#subprocess.call([sys.executable, 'C:\OSGeo4W64\bin\gdal_calc.py', '-A', 'C:\Users\jack\Desktop\New_Trash\vhFloat32.img', '--outfile=C:\Users\jack\Desktop\New_Trash\NEWCALCVH-2.img','--calc=A*(A>=1)'])


# Run gdal_calc.py to dB for HV

#subprocess.call([sys.executable, 'C:\OSGeo4W64\bin\gdal_calc.py', '-A', 'C:\Users\jack\Desktop\New_Trash\NEWCALCHV-2.img', '--outfile=C:\Users\jack\Desktop\New_Trash\HVindB.img','--calc=10*log10(power(A,2))-83'])


# Run gdal_calc.py to dB for HV

#subprocess.call([sys.executable, 'C:\OSGeo4W64\bin\gdal_calc.py', '-A', 'C:\Users\jack\Desktop\New_Trash\NEWCALCVH-2.img', '--outfile=C:\Users\jack\Desktop\New_Trash\VHindB.img','--calc=10*log10(power(A,2))-83'])


# Open Rasters using GDAL

hvRaster = gdal.Open("C:\\Users\\jack\\Desktop\\New_Trash\\hvFloat32.img")
vhRaster = gdal.Open("C:\\Users\\jack\\Desktop\\New_Trash\\vhFloat32.img")


# Get RasterBand

hvRasterBand = hvRaster.GetRasterBand(1)
vhRasterBand = vhRaster.GetRasterBand(1)


# Get Raster Mean

hvMean = hvRasterBand.GetStatistics(0,1)[2]
vhMean = vhRasterBand.GetStatistics(0,1)[2]


#### Maybe not needed
print hvMean
print vhMean

# calculate stDev

rasterList = (hvMean,vhMean)
stDev = np.std(rasterList)

#print stDev

if stDev >= 0.06:
    print "The imagery isn't acceptable"


if stDev <= 0.06:
    print "The imagery is acceptable"