Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  920 1676 (N34E27 L=267)  1169    143        3 ( 859/1631- 979/1722)  -4.8 E32-E20 N31-N36
      2  557 1655 (N24E52 L=242)   516     53        0 ( 537/1622- 577/1688) -51.3 E54-E49 N21-N26
      3  662 1645 (N25E43 L=251)   408    111        0 ( 611/1620- 711/1670)   0.0 E48-E37 N24-N26
      4  716 1545 (N20E34 L=260)   513    138        3 ( 687/1493- 740/1612) -49.4 E37-E31 N16-N24
      5 1495 1587 (N46W20 L=314)   384     74        0 (1472/1558-1527/1604)  31.0 W17-W22 N42-N47
      6  780 1486 (N18E27 L=267)   387    109        1 ( 747/1439- 807/1534) -81.9 E28-E25 N14-N21
      7  965 1389 (N17E12 L=281)   244     47        1 ( 943/1368- 992/1410)  45.0 E13-E10 N14-N17
      8  295 1278 (S04E56 L=238)   463    107        1 ( 276/1237- 317/1327) -45.0 E57-E52 S06-S01
      9  205 1210 (S10E64 L=230)  3800    212       10 ( 135/1119- 265/1298) -37.6 E77-E56 S15-S05
     10 1884 1147 (N27W64 L=358)  1074     72        1 (1845/1100-1907/1189)  71.6 W57-W69 N23-N29
     11  146 1132 (S16E72 L=222)  1806    113        5 ( 114/1091- 177/1186)  33.7 E79-E65 S18-S12
     12 1496 1098 (N13W27 L=320)   321     98        0 (1472/1051-1519/1143) -38.7 W23-W28 N10-N14
     13  594  721 (S31E21 L=273)   226     58        1 ( 571/ 699- 615/ 747) -18.4 E22-E18 S32-S29
     14 1667  569 (S12W57 L=351)   442     51        1 (1640/ 546-1693/ 592)   9.5 W54-W60 S13-S11
     15 1418  499 (S23W41 L=335)   423     61        0 (1403/ 460-1432/ 541) -45.0 W38-W43 S25-S20
     16 1531  479 (S21W51 L=345)   813    189        1 (1480/ 419-1616/ 528)  63.4 W48-W55 S26-S16
     17 1054  225 (S54W37 L=331)   538     73        0 (1034/ 201-1080/ 259) -24.0 W32-W42 S55-S50
  

    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