Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1302 1507 (N32W03 L=123)   378    110        0 (1267/1449-1330/1543) -45.0 W00-W06 N28-N33
      2 1718 1402 (N39W40 L=160)   444     73        0 (1687/1380-1751/1422)   7.1 W36-W44 N37-N40
      3 1639 1271 (N27W33 L=153)  1327    239        1 (1561/1167-1711/1333) -19.7 W25-W40 N22-N30
      4  556 1291 (S01E35 L= 85)   750    112        0 ( 488/1243- 601/1327) -12.5 E40-E31 S03-N00
      5 1506 1265 (N23W22 L=143)   448    123        0 (1472/1209-1531/1319)  82.9 W21-W23 N18-N26
      6 1784 1233 (N30W48 L=168)   376     98        0 (1774/1185-1795/1283) -80.5 W46-W48 N27-N33
      7 1398 1160 (N13W18 L=138)   302     71        0 (1358/1141-1437/1180)  45.0 W15-W19 N11-N15
      8  316 1129 (S16E49 L= 71)   293    105        0 ( 284/1094- 363/1155)  39.8 E51-E45 S18-S13
      9 1807 1018 (N18W53 L=173)  1626    210        5 (1743/ 941-1869/1085) -36.9 W45-W60 N13-N22
     10  176  923 (S31E64 L= 56)  2722    289        3 ( 130/ 815- 217/1051)  38.7 E79-E58 S38-S22
     11  791  910 (S17E10 L=110)   195     98        1 ( 773/ 861- 812/ 954) -80.5 E11-E09 S19-S13
     12 1658  785 (N00W46 L=166)   786    160        2 (1631/ 711-1680/ 863) -69.4 W43-W48 S04-N04
     13 1347  793 (S08W25 L=145)  1653    470        3 (1159/ 731-1500/ 849)  39.5 W14-W33 S15-S01
     14  188  801 (S39E65 L= 56)   669     51        2 ( 160/ 777- 211/ 826)  20.0 E71-E59 S40-S36
     15  712  695 (S32E10 L=110)   920    276        2 ( 624/ 604- 821/ 772)  57.3 E15-E04 S39-S24
     16 1712  551 (S10W63 L=184)   482     22        0 (1689/ 534-1736/ 568)  45.0 W60-W65 S11-S08
     17 1438  312 (S32W59 L=179)  2146    316        4 (1303/ 274-1568/ 343)  12.3 W45-W80 S37-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