Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1012 1835 (N59W39 L= 51)   442     67        2 ( 977/1820-1043/1853) -45.0 W34-W43 N56-N60
      2  957 1789 (N58W28 L= 40)   417     33        1 ( 935/1768- 983/1814)  68.2 W24-W31 N54-N59
      3 1572 1663 (N27W64 L= 76)   938    129        2 (1526/1632-1631/1689) -45.0 W58-W68 N23-N29
      4 1473 1510 (N22W44 L= 55)   852    170        3 (1413/1478-1539/1550) -53.1 W38-W47 N18-N26
      5 1616 1425 (N13W50 L= 62)   526    129        1 (1597/1370-1635/1490)  50.2 W46-W52 N10-N17
      6 1291 1334 (N18W24 L= 36)   477     97        0 (1245/1314-1344/1356) -45.0 W21-W27 N15-N20
      7 1103 1233 (N17W10 L= 22)   451    107        2 (1073/1196-1140/1282) -80.5 W07-W11 N14-N20
      8  553 1199 (N28E26 L=346)   517    131        2 ( 496/1167- 606/1225) -20.6 E30-E21 N25-N30
      9 1358 1049 (N00W20 L= 31)  1315    313        3 (1299/ 929-1415/1192)  40.2 W14-W27 S06-N06
     10  606 1021 (N16E25 L=346)   344    103        1 ( 581/ 978- 646/1062) -59.0 E26-E23 N13-N18
     11 1501  783 (S19W23 L= 35)   361    105        0 (1449/ 757-1548/ 802) -33.7 W19-W25 S21-S16
     12 1394  591 (S28W12 L= 24)   735    204        1 (1314/ 554-1481/ 652)  -4.4 W06-W19 S30-S26
     13 1657  594 (S37W34 L= 45)   804    104        2 (1607/ 555-1695/ 625) -18.4 W28-W37 S40-S34
     14 1544  482 (S42W22 L= 34)  1101    172        0 (1471/ 423-1618/ 532)   8.1 W15-W29 S44-S39
     15 1433  410 (S44W12 L= 23)   348     55        1 (1404/ 389-1459/ 429) -51.3 W09-W14 S45-S40
     16  506  291 (S28E67 L=305)  1873    118        2 ( 452/ 255- 547/ 326)   0.0 E77-E58 S31-S24
     17  792  282 (S34E40 L=332)   738     98        2 ( 759/ 250- 841/ 322) -81.9 E42-E37 S37-S30
  

    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