Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1276 1754 (N49W44 L=327)   535     68        0 (1246/1722-1312/1775)  14.0 W37-W49 N46-N49
      2 1405 1678 (N40W48 L=332)  1255    178        2 (1357/1619-1472/1739) -90.0 W43-W50 N34-N45
      3 1663 1599 (N28W68 L=351)   535     89        0 (1629/1562-1693/1631) -90.0 W65-W69 N24-N30
      4 1354 1601 (N36W38 L=322)   297     53        0 (1325/1573-1378/1622) -63.4 W36-W39 N33-N37
      5 1290 1549 (N34W30 L=314)   379    120        0 (1264/1506-1327/1603)  81.9 W29-W32 N30-N37
      6 1641 1485 (N22W55 L=339)   440     50        2 (1612/1464-1667/1508) -45.0 W52-W58 N20-N23
      7 1608 1339 (N14W45 L=329)   264     43        0 (1587/1310-1625/1360) -71.6 W44-W46 N12-N15
      8  850 1171 (N18E09 L=275)   849    375        1 ( 687/1111-1008/1223) -33.7 E18-W00 N11-N23
      9  364 1150 (N23E45 L=239)   678    189        1 ( 317/1061- 391/1210) -60.9 E49-E42 N17-N26
     10  234 1143 (N23E59 L=225)   357     50        0 ( 214/1121- 254/1169) -45.0 E61-E56 N21-N24
     11 1812  923 (S15W55 L=339)  5036    656        6 (1722/ 712-1925/1119)  65.4 W45-W73 S27-S03
     12  732  999 (N09E18 L=266)   405    108        2 ( 683/ 956- 776/1033) -51.3 E21-E15 N06-N11
     13 1366  890 (S07W19 L=303)   706    175        1 (1329/ 811-1413/ 956)  45.0 W15-W22 S11-S03
     14  382  463 (S21E61 L=222)   484     90        0 ( 363/ 425- 397/ 514)  30.3 E67-E55 S24-S17
     15  512  417 (S25E50 L=233)   498    104        1 ( 493/ 367- 526/ 471)  45.0 E54-E46 S29-S21
     16  579  426 (S25E44 L=240)   435     86        0 ( 543/ 390- 616/ 452)  18.4 E48-E39 S27-S23
     17 1277  298 (S49W06 L=290)   997    234        2 (1169/ 239-1357/ 356)  18.4 E05-W13 S52-S44
  

    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