EddyPro学习笔记

EddyPro学习笔记4.2版本更新内容4.1版本存在bugLicor设计了一种数据格式*.ghg这种格式里包含了需要的meta信息。如果是其他的数据格式，比如说ASCII，tob1等，需要在处理前给出meta信息。一个是快捷模式(Expressmode)，一个是高级模式(Advancedmode)。在高级模式下，需要选择各项处理采用的具体方法。可以把高级模式下设置的参数保存到express模式下，这样便于同样处理的操作。Biometdata:biologicalormeteorologicaldata生物或气象数据在project里设置timestampsreferstobeginningofaveragingperiod或者endofaveragingof需要在原始数据中查看文件显示时间到底是初始时间还是结束时间而EddyPro处理完的数据显示的时间点是结束时间点，如12:30代表的是12点到12点半的数据。在rawfiledescription设置中，把时间及records列忽略掉，选ignoreyes。如果数据不是numeric类型的，即使忽略了也要说明。Records也不是numeric类型。虽然打开时显示是10列，但是在metadata设置中应该设置成11列日期和时间算两列。因为第一列包含了日期和时间，只有半小时整时才会显示日期，其他时间显示的只有具体时间，这样的话长度不一致没法读，所以默认分成两列来读。Measurementtype仅适用于气体浓度，其他默认不需要进行设置。在metadata中要设置一下fileduration，在切割时最好都切割成30min的。对于genericbinary，需要给出一些原始数据的具体信息，例如：头文件的行数等tob1格式的数据：需要设置是ulong和IEEE4还是ulong和FP2EddyPro不能解译Ulong数据，比如时间数据通常是ulong格式，因此在rawfiledescription里面需要把时间设置成ignoreSLT(EddySoft)是EddySoft获取的特定数据格式SLT(EdiSol)是爱丁堡大学的数据获取工具获得的数据格式。Metadata数据元信息，对数据做出解释例如数据是从哪、如何获取的，用什么仪器，数据在表格中的具体位置及单位等Biometdata生物以及气象数据externalfile单个文件externaldirectory有包含气象和生物数据的文件夹但是这些数据必须按照EddyPro规定的格式组织canopyheight有效平均高度必填项displacementheight零平面位移如果不填系统会自动根据冠层高度估算Roughnesslength粗糙度如果不填系统会自动根据冠层高度估算三维超声风速仪信息中，Height必须给出，是指地面到仪器中部的距离northalignment：对于gill的三维超声风速仪，需要给出轴是否是朝北，只有gill的仪器需要给出这一项，其他不需要。Northoffset:参照当地磁北极的偏移角度，往东为正。Northwardseparation:当前三维风速仪与参考三维风速仪的距离。如果当前的在参考的北边，则为正。不用管这个项。EastwardSeparation：当前三维风速仪与参考三维风速仪的距离。如果当前的在参考的东边，则为正。不用管这个项。Verticalseparation：当前三维风速仪与参考三维风速仪的距离。如果当前的在参考的上边，则为正。不用管这个项。Longitudinalpathlength：爪子之间的距离。比如说CSAT3的距离，参照仪器的具体参数。Transversalpathlength:横向的距离。传感器的直径。Timeresponse:三维风速仪的时间响应特征，他的倒数决定了仪器能处理的大气湍流最大频率。Softwareversion：只有对Li7200和Li7500A是必须的。Height：地面到仪器中部的高度。TubelengthTubeinnerdiameterNominaltubeflowrate以上三项是针对闭路的，对于开路来说可以完全忽略。Extinctioncoefficientofwater:kw，在氪或Lyman-α湿度计中，水汽的消减系数，与参考条件下的Lambert-Beer定律的泰勒三次扩展有关。Extinctioncoefficientofoxygen:koRawFileDescription一般需要设置成ignore的列包含timestamp，linecounter等日期与时间等都不是numericmeasurementtype：仅针对气体浓度，其他变量留白或者选otherlinearscaling：一般有支持的单位就不需要再进行rescaleoutputunit：如果你进行了单位转换，那就在这里输入转换后的单位，这与下面的gainvalue和offsetvalue相对应gainvalue：斜率offsetvalue：截距Nominaltimelag：expected时滞，与三维风速仪的时滞，可根据两个仪器之间的距离估算Minimumtimelag最小时滞Maximumtimelag最大时滞左边有个选项叫rawfilesettings需要设置fileseparatorcharacter用什么隔开每个值numberofheaderrows头文件有几行这个一定要注意填写，否则会造成运行错误。BasicsettingsPageFluxaveraginginterval与numberoffilestomerge、fileduration一般平均间隔时间与fileduration相同，或者fileduration是平均时间的整数倍，不用mergeSelectitemsforfluxcomputationMastersonic:当有多个三维风速仪时，选择用哪个Fasttemperaturereading如果除了风速仪之外还有其他的高频气温观测，则可以直接用该温度替代超声虚温，这样就避免了针对超声温度的订正，比如说crosswindcorrection以及humiditycorrectionCrosswindcorrection：如果软件里没有进行此订正，在此勾选此项。一般(CSAT3)已经进行过此项校正。Flags根据某个变量的阈值把整条记录都删除掉。Advancedprocessingsettings包含四个页面处理选项、统计测试、光谱矫正、输出文件windspeedmeasurementoffsets:三维风速仪测量的风速可能有系统偏差。风速迎角订正：仅针对gill的三维风速仪，为了补偿flowdistortion效应坐标旋转：倾斜订正平面拟合使超声风速仪的坐标系与自然流线坐标系吻合，在较长的时间尺度上，例如两周或更长。可以根据不同的windsector(例如可以把360°分为12个sector，每个sector单独计算参数)计算不同的旋转角度。点击planarfittingsettings设置参数还有一个选项是planarfittingwithnovelocitybias与平面拟合类似，但是假设垂直风速的任何偏差都已经进行补偿，并强迫拟合平面过原点，也即，如果u和v等于0，则w也等于0planarfit参数设置start用于平面拟合的数据起始日期end用于平面拟合的数据结束日期这个时间段不能超过1个月，在这期间一般冠层高度与结构最好是未发生显著变化如果冠层结构与高度等未发生显著变化，2weeks的结果可以用来处理其他更长时间段的。如果规定的windsector时间越长，或者最低数据量越大，需要的时期越长。Minimumnumberofelementspersector数据太小会导致参数计算不准确。如果数据量太大可能根据就不存在此参数。Maximummeanverticalwindcomponent设置一个最大垂直风速，滤除那些过大的不实际的值Minimummeanhorizontalwindcomponent设置一个最小水平风速当水平风速过小时，迎角会受误差的较大影响，会导致最终参数不好Ifplanarfitcalculationsfailforasector如果某个sector里的参数因为任何原因，没有计算出来，选择是用临近sector的参数，或者是直接用二次坐标旋转。Setequallyspaced勾此选项会使EddyPro将360°平均分配成n个sectornorthoffsetfirstsector设置一个角度，所有的sector都顺时针从北向转此角度Turbulentfluctuations湍流波动去趋势方法选择：blockaverage简单减去平均数lineardetrending去除线性趋势做线性趋势的时间段可以不同于通量平均时间在timeconstant里设置典型的线性去趋势线性趋势时间等于通量平均时间timeconstant=0runningmean滑动平均Exponentialrunningmean类似于滑动平均但越远权重越小Timeconstant应用于lineardetrending、runningmean以及exponentialrunningmean。此参数越大，越多的低频信息被剔除。对于lineardetrending单位是分钟，对于runningmean单位是秒。Timelagscompensation：选择一种方法来补偿三维风速仪与其他高频观测之间的时滞。这种时滞一般是由于传感器之间的自然距离引起的，气流通过两个传感器时会有一定的时间差。选Constant，EddyPro会自动利用metadata中的Nominaltimelag。Covariancemaximization：利用协方差最大化计算一个窗口内最可能的时滞。窗口的大小是由最大和最小时滞定义的。Covariancemaximizationwithdefault:与covariancemaximization类似。在没有计算出最大协方差时，利用Nominaltimelag。Automatictimelagoptimization：选此选项，需要在Timelagoptimizationsetting中进行设置，他会计算出Nominaltimelags以及Plausibilitywindows，并在原始数据处理中进行应用。TimelagoptimizationsettingsTimelagfilenotavailableStartEnd起始和结束时间，时间段一般在1到2个月。Plausibilityrangearoundmedianvalue这个范围是由中位数±n个MAD（MAD是absolutedeviation，绝对偏差）Watervaportimelagasafunctionofrelativehumidity水分的时滞把RH分级，例如分为10个等级，则0-10%为等级1。。。如果选择1个等级，则会把水分与其他气体一样处理，对于开路或者管子非常短的闭路是合适的。惰性气体PassivegasesMinimum(absolute)CO2flux低于此CO2通量值对应的CO2时滞在计算timelag时将不被考虑，这样选择足够大的通量值，对应较好的湍流交换情况。TimelagsearchingwindowsMinimumMaximum设置最低和最高时滞，启动优化程序。由最低和最高时滞决定的这个范围一定要足够大，使其可能涵盖所有可能的时滞。密度效应补偿Compensationfordensityfluctuations(WPLterms)如果是混合比或者有所需资料可以将原始资料直接转化成混合比的话，就不用进行密度效应校正了。Burba校正如果是7500的话，加上仪器加热效应产生的感热通量部分。有两种，一种直接根据温度计算，另一种是多元回归，根据温度，短波入射辐射，长波辐射以及风速。快速傅里叶变换Fastfouriertransform对时间序列数据进行频域分析，