Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1645 1592 (N18W61 L=174)  2372    297        4 (1600/1474-1720/1673)  14.5 W52-W80 N12-N22
      2 1203 1616 (N33W30 L=143)   309     95        0 (1162/1599-1250/1633) -18.4 W26-W32 N30-N34
      3 1204 1514 (N27W25 L=138)   645    167        2 (1136/1463-1277/1547)   4.8 W18-W30 N25-N27
      4 1239 1397 (N19W22 L=135)   270    126        0 (1179/1361-1292/1421) -50.2 W19-W24 N15-N21
      5  717 1219 (N22E13 L=100)   907    241        2 ( 664/1137- 803/1298) -67.4 E15-E09 N15-N27
      6 1877 1148 (S12W62 L=175)  1954    261        3 (1846/1045-1914/1271)  42.5 W56-W70 S17-S05
      7 1629 1102 (S07W38 L=151)  1686    394        3 (1569/ 914-1678/1254)  66.0 W30-W42 S16-N02
      8  340 1221 (N31E42 L= 72)   932    176        0 ( 274/1190- 415/1249)  -4.1 E48-E34 N28-N32
      9  543 1150 (N22E26 L= 87)   910    253        2 ( 437/1077- 643/1196) -37.6 E33-E20 N16-N26
     10  858 1129 (N14E07 L=106)   177     76        0 ( 817/1115- 885/1149)   0.0 E09-E04 N12-N14
     11  351 1039 (N20E43 L= 71)   410    133        1 ( 286/1015- 413/1069) -28.6 E48-E37 N17-N23
     12 1898 1016 (S20W62 L=175)   728     94        2 (1876/ 971-1917/1068)  45.0 W58-W65 S22-S16
     13 1743  820 (S27W42 L=155)  1827    256        4 (1626/ 751-1845/ 887)   2.6 W31-W53 S28-S23
     14 1445  800 (S19W19 L=132)  1018    235        3 (1341/ 744-1545/ 867)   0.0 W11-W26 S20-S17
     15  738  776 (S03E22 L= 91)   693    310        2 ( 624/ 698- 860/ 859) -53.1 E27-E16 S09-N03
     16  221  627 (S01E69 L= 45)  1033     96        3 ( 196/ 589- 247/ 679) -78.7 E72-E65 S03-N02
     17  307  553 (S07E62 L= 51)  5110    277       12 ( 228/ 429- 378/ 665)  54.5 E74-E53 S14-N00
     18  436  531 (S11E51 L= 62)   979    168        4 ( 403/ 461- 465/ 598)  71.6 E54-E48 S15-S06
  

    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