Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  921 1795 (N52E26 L=193)   336     75        2 ( 898/1767- 939/1815) -39.8 E29-E23 N49-N54
      2  427 1437 (N18E48 L=171)   801    209        1 ( 409/1344- 448/1518) -70.0 E51-E45 N12-N23
      3 1322 1456 (N32W15 L=234)   288     76        2 (1295/1424-1358/1478)  53.1 W12-W17 N29-N33
      4 1153 1396 (N25W03 L=223)   363    137        1 (1134/1334-1174/1460)  90.0 W02-W04 N21-N29
      5 1428 1414 (N30W23 L=242)   472     91        1 (1386/1371-1468/1449) -26.6 W19-W25 N27-N32
      6  256 1170 (S00E57 L=162)   404     74        2 ( 234/1140- 279/1202)  90.0 E59-E54 S02-N01
      7 1756  868 (N00W53 L=272)  5192    865        9 (1679/ 559-1832/1186) -78.4 W46-W64 S19-N20
      8 1695 1114 (N14W44 L=264)   554    111        2 (1654/1077-1734/1143) -74.1 W41-W48 N12-N16
      9  247  934 (S15E55 L=164)  1066    152        3 ( 221/ 861- 277/1002) -60.9 E57-E51 S19-S10
     10  212  830 (S22E59 L=161)  2286    343        7 ( 165/ 708- 279/ 912) -56.3 E64-E51 S29-S16
     11  675  764 (S20E19 L=200)   338    124        0 ( 623/ 725- 733/ 807) -24.0 E23-E14 S21-S17
     12 1452  672 (S16W33 L=252)  2150    544        7 (1333/ 498-1607/ 799)  72.3 W27-W43 S28-S06
     13 1291  779 (S11W20 L=239)   290     88        1 (1239/ 762-1332/ 794)   9.5 W16-W22 S12-S09
     14 1521  655 (S16W38 L=257)   672    172        1 (1466/ 612-1587/ 686)  40.6 W34-W43 S19-S13
     15  486  573 (S36E34 L=185)   430    119        2 ( 444/ 548- 525/ 613)  32.0 E38-E30 S37-S32
     16 1505  420 (S32W46 L=265)   479     99        1 (1464/ 403-1546/ 448)  26.6 W43-W49 S33-S29
  

    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