常見探地雷達數(shù)據(jù)格式

GSSl?DZTHerearethetwoshucihiresusedintheDZTfileheader.unsignedsec25/*ssccnd/2unsignedmiunsignedsec25/*ssccnd/2unsignedmin6/*minuteunsignedhour5/*hourunsignedday5/*dayunsignedmonth4/*mnthunsignedyear7/*year-1980structDztDateStructbits)(0-29)(0-59)(0-23)(1-31)(1-12；*/*/*/*/[二」anr2=Feb.etc?(0-127j1900-2107)*/iftightlypacked*//*N:?te:thefbllowingstructureincorpiratesthechangesnotedinthe' preIiminarydraftforextendedRADANfileh&aderformatas* faxe?dto」曰FfLuciusbyLeoGaIinovsky(GSSI)1/11/94.strLEtDztHdrStruct{ unsignedshortrh_tag;/*OxONff,wliereN=rh_nchan-1(0-15)*/unsignedshortrh_dataj/*offsettodata(1024*rh_nch3n)*/unsignedshortrh_nssmp;/*samplesperscan(2-65S35)*/unsignedshortrh_t?itsj/*bitsperdatawc?rd(S.16,32,64)*/shortrh_zerojZ*binaryoffset(-128.-327G8retcJ?/floatrh_sps;/*scanspersecond*/floatrh_spm;/*scanspermeter*/floatrh_mpni:/*meterspermark*/floatrh_position;/*position(ns)*/floatrh_range：Z*range(ns)*/unsignedshortrh_np?asBi/*scansperpassfor2Dfiles*/structDztDateStructrh_createj/*datecreated*/structDztDateStructrh_jncdif：/*datemodified*/unsignedshortrh_rgain;/*offsettorangegainfbnction*/unsignedshortrh_nrgainj/*sizeofrangegainfunction*/unsignedshortrh_textj/*offsettotext*/unsignedshortrh_ntexti/*sizeoftext*/unsignedshortrh_procj/*offsettoprocessinghistory*/unsignedshortrh_nproc：/*sizeofprocessinghistory*/unsignedshortrh_nchan;/*nunberofchannels*/floatrh_epsri/*averagedielectricconstant*/floatrh_top;/*toppDSitioninmeters*/flostrh_depth;/*rangeinmeters*/charreserved[31]charrh_dtype：/*bits765432100unsigneddata1signeddata08-bitint116-bitint132-bitint164七itint132-bitfloat164-bitfloat(USGSextna)

XXX notused/*rh_dtypadoesnotappearinversionsbefore5：2*/char rh_antnam^[14]:/?antennanaiw($g.3105(SOOMHz))和urb^ignedshort rh_e：h^nni3^k■ /?setivsehsnIirs-skformatib0k530Xf通h宕「臼Kh-asfolIcwingbitsset:bit§!7￡54321a11chanriel112channels1113chann@la11114chanrtels000Qnotusedchar rh_name[12]:/?thisfilenamewithc-utDOSextensionunsignedshort rh_chksum;/*checksumforheader:NOTE:saeprcgrambelowforproperwaytocalculateheaderchecksum(essentially,zerooutoldchecksumvalue(ccpyitfirst)and<iddupvalues,inthe-hG^deras.ifit縫宵「?&512unsignedshorts,thenplaceGhwk^umvalueintorli_Ghk-sum??/八128byte-sto2“4/char variable[89&]; /?rangegain,corrments,andprocessinghistory"}; /?1024bytesiftightlypacked?/WARNING!InolderDZTfiles,therl\_datafieldmaynotactuallybetheoffsettodatainmuIti-chann^Ifiles!Itisprotoblyb^sttoreadinaIItheheadersichleavesthefilepninterat:thestartoFthefirstscanorreadintheFirstheadertodiscoverthe*numberofchannelsthemrewindtheFile*andmovethefilepointer"numberc-fchanneIstimessizwfleaderstructure)**?HereisaFunctionthatprintsoutthecontentsofaDZ.Tfileheader.11aIsoshwshwthechecksumiscalculatedaridb>wtoreadtherangegain?? prjntOneDzthte^cfer()Printouttheccnt^ntsofoneDZTfil&headertoonc^itF^itdevice(stdout,9td?r「.adiskFiIp,etc?)??Parameters:intheaderjDiytes-numberofacXuaIbytesinfileheader,512or1024intnum_traces ?numt:$「ofradar-scans,<=0ifunkncwicharFuj和lenamm-"stdout"stderr11,oraDOSfiIchar*dzt_fi ?DOSnan>=?ofthe?DZTfileheader(%) longtostructDztHdrStruct*hdrPtrl-addressesoftheheaderffequires:<ccnio.h>.<stdio.h>r<string?h>,"gpr—io?h"?C<i11s:fopsn.fprintf,getchi,strstr^ConvertProcHist2.R?turn^i0onguomsg..>0onerror(offsetintoanarr<iyofstring^-)?*?Author:JeFfLuciusUSGSBranchofCeophysics GoIden,8?Data:D9CG<nbGr*1,1Q95■/intPrintOn^DztH*?ader(intheaderJongnLfn-trac9S.,charfout_filenani?,cha廠4dzt_.fiIename,struct:DztHdrStruct*bdrPtr){charant_freq[10].n^A'Str[32]:unsignedshortrgjDreaks=0;unsignedshortchecksum:unsignedshortustemp;inti.「jb廠Ita;FILEStream=NULL;staticexternextern?xternicharbig_buff[1O24]staticexternextern?xternicharbig_buff[1O24]:/?removesFromstackandplacesindataeonGtcharconstcharconstchar*/?ant.number[]: /? in*ant_nan>?[]: /* infm5nth_abbr[] ： /? in??dzt_io.hM*/-,dzt_iQ.hM*/?dzt_io.hMAssignorop?nstreamAssignorop?nstream?/if (strstr(out.Filename,"stdout"))stream=stdout;el if(strstr(out_filenam目,"stde?廠「“5jstr^am=stde-rr;與I3 stream=fopen(c<it—fiIename,"wt*');if(stream==NULL)return1:/*Printtheheecterinformation?/?fprintf(stream,''SIR-iOAheaderforFile%s:\n",dzt_-fiIenarr??);fprintf(str^m.-,numt*?rofheaders=%d、rT“hd「Ft「?^h-nchan);fprintf(stream,e,Fileheadersize=td\n"^header_bytes):if(nunL.traces>0)fprintFfstreani,"numberofGPRtraces=Hd\nB,nuni_traces):fprintf(str-^m,"SIR-10AheaderID=ftxX04x\n“.hd「Pt「?A「h_tQg》；fprintf(str^m,"offsettodat-a =lhd\n-,fhdrPtr->rh_datA):

fprintf(strearn."sampIespersean =4hu\n，!HhdrPtr->rh_nsamp);fprintf(strearn#Hbitsperdateword=4hu\n，!HhdrPtr->rh_bits)jfprintf(stream,!fbinaryoffset =%hd(0kXO4X)\n!!,hdrPtr->rh_zeroP(unsignedshort)hd廠Ptr->rh_zero);fprintf(strearnH”scanspsrsecond=HfXn"HhdrPtr->rh_sps);fprintf(strearn『fprintf(strearn『"scanspBrmeterfprintf(strearn『!，niterspermarkfprintf(strearn*!，p:?sition(ns)fprintf(strearn"range(ns)fprintf(strearn.!，scanspe「passif(hdrPtr->rh_create.mnth=0)fprintf(stream,pcreatedoteeIse=%fXn"..hdrPtr->rh_spoi);=ftf\n"..hdrPtr->rh_mp<n);=l.fXn"..hdrPtr->rh_position)i=IsfXn：.hdrPtr->rh_rangs);=%hu\n",hdrPtr->rh_npsiss);=nodsteXn");fprintf(stream丿scr&atedate =feXd,%d(%d:^02d:^O2d)\n!!.mnth_8bbr[hdrPtr->rh_create.month].hdrPtr->rh_create.day”hdrPtr->rh_create.year+19B0?hdrPtr->rh_create.hour’hdrPtr->rh_crs8te.min,hdrPtr->rh_crsate、sec2*2)jf(hdrPtr->rhjiedif.month==0)fprintf(9treamBmodificationdate=nodate\nu);eIsefprintffstream?Bmodificationdate=fe%dfi%d(fel:XO2d?W2d)\n!!.n>Dnth_abbr[hd廠Ptr->rh_mcdif.in^nth].hdrPtr->rh_nKidif.day,hdrPtr->rhjn>Ddif.year+1980*hdrPtr->rhjn>Ddif.hourrhdrPtr->rhjiedif.min.IhuXn",lidrPtr->rh_rgain);thu\nuIhuXn",lidrPtr->rh_rgain);thu\nu,hdrPtr->rh_nrg8in);%hu\nu,hdrPtr->rh_teKt);UhuVn**.hdrPt廠?>rh_ntext);￡hu\n'',hdrPtr->rh_proc);4hu\nu,hdrPtr->rh_nproc);fprintf(str&am#Hoffsettorangegainfprintf(strearn/psizeofrangegsinfprintf(strearn totextfprintf(strearnHnBizeoftextfprintf(strearn."offsettoproc.hist,fprintf(strearn#nBizeofprocahistsif(Btream==stdout11stream==stderr){fprintfCstream,"Pressakeytocontinue...getch0;/*fprintf(stream..B\rputeC111};fprintf(strearn.!，numtefprintf(strearn.!，numte「ofchanneIsfprintf(strearn.!，sve,dieI.constantfprintf(strearn『"toppos.fprintf(strearn"rangeinfprintf(stream^"rsserv&dant_freq[0]=D;for(i=0j日nt_number[i]jinmetersnitersi++)=thu\n!!,hdrPtr->rh_nchan);=ftfXn"..lidrPtr->rh_epsr);=%fXn"..hdrPtr->rh_top);=l.fXn"..hdrPtr->rh_depth);=te\n"..hdrPtr->reserved);{if(&trstr(hdrPtr->rh_antna隕{if(&trstr(hdrPtr->rh_antna隕break: /*outofforIoop*/}fprintf(strearn土"antennaname =fc"..hdrPtr->rh_antnan>9)jif〔日nt_freq[0])fprintf(stream「(%s)"‘a(chǎn)nt_freq)\fprintf(&trearn,"X.n"):fprintf(str&am#"activechannelsmask=0x%Q4X\n".hdrPtr->rh_chanmask);fprintf(&trearn,"thisfiIenan>9 =teXn"hdrPtr->rh_name);/*Cbleulatechecksumforheader*/if(hdrPtr->rh_chksumI=0){fprintf(streamHChecksumforh&acfer=%liu(0x^D4X)\n'HhdrPtr->rh_chksum,hdrPtr->rh_chksum);ustemp=hdrPtr->rh_chksumjhdrPtr->rli_chksum=Ojchecksum=0jfor(i=0-i<S12ii++)checksijm+=*((junsignedshort*)hdrPtr+ijjfprinff(strsam?pcaleu checksum=%liu(OxitMX)[headercorrupted]\n"checksum,checksum,,(checksum==ustemp)?"isnot":"is");hdrPtr->rh_chksum=ustemp;}if(hdrPtr->rli_rgain!=0&&hdrPtr->rh_nrgain!=0){rg_breaks=?(unsignedshort?)((cliar4)hdrPtr?hdrPtr->rh_rgain):fprintf(stream,*numberofrgbreaks=$Uid\n-,rgj>reaks):fprintf(stream,"rangegain:scansampledb\n"):rg_break_.de!ta=(hdrPtr->rh_nsamp-1)/(rgj>reaks-1):for(i=0;i<rgj>reaks;i{fprintf(streamrH %3d 篤7.3Cn”.rqj^reak.deIta*ir?(flc?at?)((char4)hdrPtr?hdrPtr->rh_rgain?2?4*i));}}if(streams=stdout||stream==stderr){ fprintf(stream,*Pressakeytocontinue???");getch();/?fprintf(stream,*\r '「“)；?/puts(““)；}if(hdrPtr->rli_text)fprintf(stream,?coiwnents=\n%s\n°,(char*)hdrPtr?hdrPt「?>「h_text):elsefprintf(stream,"nooomnantsXn");if(hdrPt「?>「l\_nproc)hbt?:ThesourcecedeforConvertProcHist2()isinthefileDZT_IO.C{ ConvertProcHist2((int)sizeof(big_buff),bigj>uff,(int)hdrPtr->rh_nproc,((char?)hd「Pt「?hdrPtr->rh_proc)):fprintf(streamz*ProcessingHistory:\nM):if(strlen(bigj>uff))fprintf(stream, biqj^uff):else fprintf(stream,M\tnone\n\nM):}elsefprintf(stream,?noprocessinghistory\nM):if (strstr(out_fiIename,Mstdout"));elseif(strstr(out_filename,"stderr"));else fclose(stream):return0;}Sensors&Software一DT1/HD/?Thisstructureispresentteforeeverytraceinthe4.011datafile?/structSsTraceHdrStrict?traceidentificationheader?/flcattrace_jium:/*tracesequencenumterwithinIine?/flcatposition;/?positionalongtraverse?/floatnum.samples;/*samplespertrace?/floateIevation:/?elevaticndataifavaiIable?/floatunass5:/*unassigned*/floatnumj>ytes:/*bytespersample;always2forRev?3firntA'are?/floataux.tracejium:/*tracenumteragain?/floatnum_stacks;/*numberofstacksusedtogettrace?/floattirr>?_windw:/?timewindsinns?/floatunasslO;/*unassigned?/floatunass!1:/*unassigned?/floatunass!2:/*unassigned?/floatunass13;/*unassigned?/flcatunass!4:/?unassigned?/flcatunass15;/?unassigned?/flcatunassIG;/?unassigned?/floatunass!7:/*unassigned?/floatunass18:/*unassigned?/floatunass19;/?unassigned?/floatunass20:/?unassigned?/floattin^_zero_adjust:/?wheresamplex=samplex?adjustment?/floatzero_flag;/?0=dataOK,1=zerodata?/

floatunass23;floattin>?:floatconment.flag;charcomn>?nt[28]floatunass23;floattin>?:floatconment.flag;charcomn>?nt[28]:/*secondspastmidnight?//*1=comn^ntattached?//*optimalcharactercoiwnents?//*128bytesiftightlypacked?/SEG一SEGYorEBCDIC-cededblock?//*theSEG-Yreelidentificationheader?/orEBCDIC-cededblock?/charcomn>?nt[3200];/?3200-bytetextarea?/longintjcbid;/*jobidentificationnunber?/longintIino;/?Iinenumber(onlyoneIineperr?l)?/longintreno:/?reelnumber?/shortintntrpr;/?numterofdatatracesperrecord*/shortintnart:/?numterofauxiIiarytracesperrecord?/shortinthdt:/?sampleintervalinmicrosecsforthisrwl?/shortintdto:/?samefororiginalfieldrecording?/shortinthns;/?numterofsarrplespertraceforthisreel?/shortintnso;/?numterofsarrplespertracefororiginalfieldshortintformat:/?datasampleformatcede:=fIcotingpoint(4bytes)=fixedpoint(4bytes)=fixedpoint(2bytes)=fixedpointw/gaincode(4bytes)?/shortintfold:/?CDPfoldexpectedperCDPensenble?/shortinttsort;/?tracesortingcode:=asrecorded(nosorting)=CDPensenble=singlefoldcontinuousprofile=horizontaIIystacked?/shortintvseode:/?verticalsumcode:=nosun=twosum???N=Nsum(N=32,767)?/shortinthsfs:/?swepfrequencyatstart?/shortinthsfe:/?frequencyater>J?/shortinthslen;/?length(ms)?/shortinthstyp:/?swemptypecede:=Iinear=parabolic=exponentiaI=other?/shortintschn:/?tracenumterofswempchanneI?/shortinthstas;/?sweptracetaj?rlength(msec)atstartiftapered(thetaperstartsatzerotimeandiseffeetiveforthislength)?/shortinthstae;/?sweptracetap^rlength(msec)atend(theendingtaperstartsatsweplengthminusthetaperlengthatend)?/shortinthtatyp:/?traceta|?rtyp^code:=Iinear=cos?squa「md=other?/shortinthcorr;/?correlateddatatracescode:=no=yes?/shortintbgrev:/?binarygainrecoveredcode:=yes=no?/shortintrevm;/?ampIituderecoverymethedcode:=none=sphericaIdivergence=AGO=other?/shortintmfeet:/?n>?asur€<nentsystemcode:recording?/=meters=feet?//?binary-codedblock?/1337.5to1337.5to22.5degrees222.5to67.5degrees367.5to112.5degrees4112.5to157.5degrees5157.5to202.5degrees6202.5to247.5degrees7247.5to292.5degrees8293.5to337.5degrees?/shortinthunass[170]: /?unassigned"V/*3600bytesiftightIypacked?/shortintvpoI:codeseismicsignallagspilotby=increasegeophonene^ative=increasegwpb^nepositive/?vibratorypolarityinpressureorupwardcasemovemgntgivesnunberontapeinpressureorupwardcasemovementgivesnunberontape?/code:/*theSEG-Ytraceidentificationheader?/stri>3tSegyTraceHdrStructlonginttrad:/?tracesequencenumberwithinlire?longinttracr;/?tracesequencenumberwithinreel?longintfldr;/?fieldrecordnunber?/longinttracf;/?tracenumt-?rwithinfieldrecord?/longintQP；/?energysourcep^intnumter?/longintcdp:/?CDPensemblenunber?/longintcdpt:/?tracenumt-?rwithinCDPensemble?/shortinttrid:/?traceidentificationcode:=seismicdata=dead=dunmy=tirr>?break=uphole=sw*p=timing=aviterbreak—,N=optionaIuse(N=32,767)FollowingareCAPidflags:=autocorrelation=Fcuriertransformed-nopackingxr[0],xi[0] xr[N-1],xi[N-1]=Fouriertransforn^d-unpackedNyquistxr[0]zxi[0] xr[N/2],xi[N/2]=Fc<iriertransformed-packedNyquistevenN:xr[0],xr[N/2]rxr[1],xi[1] xr[N/2-1],xi[N/2-1](notetheexceptionaIsecondentry)cddN:xr[0]/xr[(N.1)/2]zxr[1],xi[1] xr[(N-1)/2-1],xi[(N-1)/2-1]rxi[(N-1)/2](notetheexceptionalseccod&lastentries)=Cc<nplexsignalinthetirr>?domainxr[0]zxi[0] xr[N-1]zxi[N-1]=Fcuriertransformed-ampIitucfe/phasea[0].p[0] a[N.1],p[N-1]=Cc<nplextimesignal-ampIitucfe/phasea[0].p[0] a[N.1],p[N-1]=RealpartofcompIextracefrc<n0toNyquist=ImagpartofcompIextracefrc<n0toNyquist=ArrpIitudeofcompIextracefrc<n0toNyquist19=Phaseofcomplextracefrom0toNyquistshortintnvs;=Wavenumbertimedomain(k-t)=Wavenumberfrequency(k-ome^a)30=Depth-Rang?(z-x)traces=Seismicdatapackedtobytes(tysupack1)=Seismicdatapackedto2bytes(bysupack2)200=GPRdata?//*numterofverticallysunmedtraces(seevseedeshortintnhs:/?inreelheaderstructure)?/numt-?rofhorizontaIIysumrr>?dtraces(seevseodeshortintduse:/?inreelheaderstructure)?/datause:longintoffset:/?=prcducticn=test?/distancefromsourcepointtoreceiverlongintgelev;/?group(negativeifoppositetodirectioninwhichtheIinewasshot)?/receivergroupeIevationfromsealevellongintselev;/?(abovesealevelispositive)?/sourceelevationfromsealevellongintsdepth：/?(abovesealevelispositive)?/sourcedepthbeIa*surface(positive)?/longintgdel:/?datumelevationatreceivergroup?/longintsdel:/?datumelevationatsource?/longints^dep:/?waterdepthatsource?/longintg^dep:/?waterdepthatreceivergroup?/shortintseaIel;/?scalefactorforprevious7entriesshortintsea18；/?withvalueplusorminus10tothep*E「0,1,2>3,or4(ifpositive,multiply,ifnegativedivide)?/scalefactorforr?xt4entrieslongintsx;/?withvalueplusorminus10tothe0,1z2>3,or4(ifpositive,multiplyzifnegativedivide)?/Xsourcecoordinate?/longintsy;/?Ysourcecoordinate?/longintgx;/?Xgroupcoordinate?/longintgy;/?Ygroupcoordinate?/shortintcounit;/?coordinateunitscode:shortintweveI;/?forprevicusfourentries=Iength(netersorfwt)=secondsofarc(inthiscase,theXvaluesarelongitud?andtheYvaluesarelatitude,apositivevaluedesignatesthenumberofsecondseastofGreenwichornorthoftheequator?/gatheringvelocity?/shortintswevel:/?subA^atheringvelocity?/shortintsut;/?upholetirr>?atsource?/shortintgut：/?upholetirr>?atreceivergroup?/shortintsstat;/?sourcestaticcorrection?/shortintgstat；/?groupstaticcorrection?/shortinttstat;/?totalstaticapplied?/shortintlaga:/?lagtirr>?Aztin>?inmsZtwewnendof240bytetraceidentificationheaderandtimebreak,positiveiftimebreakoccursafterendofheader,tin>?breakisdefinedastheinitiaticnpulse*hichmayberecordedonanauxiIiarytraceorasotherwisespecifiedbytZrecordingsystem?/shortintlagb: /?lagtirr>?B,tin>?inmsZtweenthetiir>?breakandtheinitiationtimeoftheenergysource,maybepositiveornegative?/shortintdelrt;/?delayrecordingtime,timeinmsbet^wninitiationtimeofer^rgysourceandtimewhenrecordingofdatasamplesbegins(ford*pAbaterwrkifrecordingdoesnotstartatzerotime)?/$hortirrtmutG.：shortinxmute-;unsignadshortintunsignedshortintshortintgain;/? 11m^--st.art?//*mut?time--&n<l*/ns:/?numberofsamplesinthistrace?/dt;/?sampIeinterval;inmicseconds?//?gaintype*offieldinstruniGiitGcede?:1=fixedshortintigc;shortintigi;shortintcorr;二binary=flcatingpointN=optionaltee*/嚴(yán)imstrun>9ntgainconstant$hortirrtmutG.：shortinxmute-;unsignadshortintunsignedshortintshortintgain;/? 11m^--st.art?//*mut?time--&n<l*/ns:/?numberofsamplesinthistrace?/dt;/?sampIeinterval;inmicseconds?//?gaintype*offieldinstruniGiitGcede?:1=fixedshortintigc;shortintigi;shortintcorr;二binary=flcatingpointN=optionaltee*/嚴(yán)imstrun>9ntgainconstant?//*instrunentearlyorinitialgain?/八correlated:1=noshortintsfs;GhortirrtGfe:shortintslen:shortintstyp:2=?//*sw?旨pfrequencyat^tart?//?mwxpfrequencyatenJ?//*sweeplengthinms?/八svweptypecede:=Iinear=parabolic=exponential=other*/shortintshortintshortintshortintGhortintshortinxGhortintshortintshortintshortintshortintshortintshortintshortintshortintshortint5hortint fsepfG—I:;巾^1f>?ff>