Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1417 1642 (N45W35 L=131)  1413    207        2 (1374/1571-1459/1693)  42.0 W29-W41 N39-N48
      2  940 1477 (N32E06 L= 91)   340     94        1 ( 925/1438- 961/1522) -71.6 E07-E04 N29-N35
      3 1731 1463 (N30W57 L=153)   453     89        1 (1705/1440-1755/1481) -26.6 W54-W59 N28-N31
      4 1559 1317 (N21W37 L=134)   533    128        0 (1506/1269-1617/1366) -40.6 W33-W41 N18-N24
      5 1562 1226 (N15W36 L=133)   591    178        1 (1493/1202-1628/1261)   5.7 W31-W41 N14-N17
      6 1112 1166 (N12W05 L=102)   373     91        0 (1092/1120-1137/1220) -71.6 W04-W07 N09-N15
      7  604 1131 (N10E27 L= 70)   252     71        2 ( 569/1108- 640/1153) -21.8 E29-E24 N08-N10
      8 1574 1096 (N07W36 L=133)   411    119        0 (1517/1073-1624/1128) -20.6 W32-W40 N05-N09
      9  507  996 (N01E34 L= 63)   417    109        0 ( 456/ 969- 562/1028)  14.0 E37-E29 N00-N03
     10  399  701 (S18E44 L= 53)   339     64        1 ( 385/ 664- 412/ 736)  78.7 E45-E42 S20-S15
     11  727  625 (S22E20 L= 77)   668    142        1 ( 668/ 599- 801/ 662) -12.5 E23-E14 S23-S19
     12 1587  622 (S23W40 L=137)   549     79        2 (1556/ 587-1616/ 659)  68.2 W38-W42 S25-S20
     13  491  591 (S25E38 L= 59)   546    132        0 ( 442/ 573- 551/ 611)   6.3 E42-E33 S26-S23
     14 1126  408 (S38W08 L=105)   438    107        3 (1105/ 369-1145/ 449) -81.9 W06-W09 S41-S34
     15  728  376 (S41E24 L= 73)   444     78        0 ( 699/ 354- 750/ 410)  38.7 E27-E22 S42-S38
     16  995  342 (S43E02 L= 95)   875    124        1 ( 938/ 310-1038/ 387)   0.0 E07-W01 S46-S39
  

    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