Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1146 1700 (N47E03 L=118)  1053    269        3 (1108/1621-1202/1786)  90.0 E05-W00 N41-N56
      2  622 1721 (N39E46 L= 75)   529    161        0 ( 586/1679- 646/1749) -24.0 E51-E42 N36-N40
      3  381 1658 (N30E64 L= 56)  1020    164        2 ( 356/1628- 401/1697) -51.3 E66-E61 N28-N33
      4  541 1542 (N25E43 L= 78)  1857    289        5 ( 413/1471- 653/1665)  51.3 E52-E38 N19-N35
      5 1515 1583 (N45W31 L=152)   523     84        1 (1487/1552-1542/1614)  45.0 W27-W34 N41-N47
      6  578 1408 (N17E35 L= 85)   380    143        1 ( 517/1371- 619/1461)  39.8 E39-E33 N14-N20
      7 1489 1424 (N32W27 L=148)   309     74        2 (1473/1385-1504/1457) -78.7 W25-W28 N29-N34
      8  714 1249 (N09E22 L= 98)   291     84        1 ( 694/1215- 742/1282) -53.1 E23-E20 N07-N11
      9 1204 1215 (N14W08 L=129)   485    124        1 (1157/1176-1243/1271)  80.5 W05-W10 N11-N17
     10 1134 1158 (N10W05 L=126)   235     57        1 (1103/1134-1160/1177) -14.0 W02-W06 N08-N10
     11  243  955 (S15E55 L= 66)   752    239        0 ( 209/ 858- 262/1066)  77.0 E58-E52 S21-S08
     12 1304 1011 (N03W17 L=138)   229     80        0 (1289/ 963-1320/1055) -68.2 W15-W18 N00-N05
     13 1653  907 (N02W43 L=164)   434     79        2 (1632/ 868-1673/ 947) -53.1 W41-W45 N00-N04
     14  833  887 (S11E10 L=111)   335    156        1 ( 746/ 870- 905/ 908)  16.7 E15-E05 S11-S08
     15 1205  789 (S11W14 L=135)   697    193        2 (1161/ 715-1248/ 888) -63.4 W11-W17 S15-S05
     16 1630  731 (S09W45 L=166)   362    119        1 (1602/ 662-1662/ 783)  82.9 W43-W46 S13-S05
     17 1739  692 (S10W57 L=177)   406     68        1 (1714/ 662-1757/ 731)  76.0 W54-W58 S11-S07
  

    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