Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  714 1680 (N36E38 L= 65)   483     87        2 ( 681/1653- 747/1720) -24.0 E43-E34 N34-N38
      2 1402 1257 (N20W20 L=124)   283     70        1 (1366/1231-1436/1285) -21.8 W17-W22 N18-N21
      3 1851 1217 (N25W62 L=166)  2301    191        7 (1764/1181-1936/1261)   0.0 W51-W83 N22-N26
      4 1633  939 (N04W41 L=144)   274     49        1 (1606/ 919-1656/ 963)  14.0 W38-W42 N02-N05
      5 1210  647 (S20W18 L=122)   347    151        0 (1194/ 569-1233/ 712) -73.3 W17-W20 S25-S15
      6  955  496 (S35W05 L=108)   701    143        1 ( 929/ 432- 982/ 566) -53.1 W01-W08 S39-S29
      7 1211  513 (S29W22 L=126)   259     45        1 (1183/ 493-1233/ 536)  63.4 W20-W23 S30-S26
      8 1299  471 (S30W30 L=134)   322     79        0 (1257/ 447-1329/ 501)  53.1 W27-W32 S32-S27
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y