Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  955 1641 (N32E13 L= 49)   290     72        1 ( 915/1618- 988/1659) -18.4 E16-E10 N30-N32
      2  284 1487 (N16E61 L=  1)  5454    347       11 ( 199/1339- 384/1610)  74.5 E76-E53 N06-N24
      3  361 1482 (N16E52 L= 10)  1660    189        5 ( 303/1403- 422/1553)  63.4 E56-E48 N10-N20
      4  176 1405 (N11E72 L=350)   900     64        3 ( 156/1364- 198/1441) -78.7 E76-E69 N07-N12
      5 1868 1371 (N29W63 L=126)   875    136        0 (1829/1318-1902/1431) -36.9 W59-W68 N25-N32
      6  417 1232 (N01E41 L= 21)   385     84        1 ( 371/1201- 458/1253)   0.0 E44-E37 S01-N01
      7 1697 1156 (N11W42 L=104)   770    193        2 (1663/1066-1736/1248)  71.6 W39-W45 N05-N17
      8 1842 1143 (N13W55 L=118)  1141    136        2 (1797/1081-1874/1206) -45.0 W50-W59 N09-N16
      9  492  839 (S22E33 L= 30)  2597    694        3 ( 313/ 641- 675/1146)  49.7 E51-E21 S35-S02
     10  164 1061 (S11E63 L=359)  1753    177        7 ( 102/1018- 219/1110)  10.6 E72-E56 S13-S07
     11  374 1003 (S13E42 L= 20)   419    131        1 ( 332/ 971- 414/1028) -15.9 E45-E38 S15-S11
     12  237  967 (S16E55 L=  8)   372     52        1 ( 217/ 943- 258/ 998)  71.6 E56-E52 S17-S14
     13  742  930 (S14E16 L= 46)   298    120        2 ( 704/ 885- 770/ 982) -45.0 E18-E13 S16-S11
     14 1379  904 (S08W23 L= 85)   422    154        0 (1344/ 835-1405/ 974)  83.7 W21-W23 S12-S03
     15  210  875 (S22E59 L=  4)  1974    349        2 (  84/ 832- 371/ 916)   7.0 E83-E42 S23-S18
     16 1885  836 (S03W66 L=128)   547     69        0 (1859/ 799-1905/ 862) -15.9 W61-W69 S05-S02
     17 1341  730 (S19W24 L= 86)   847    324        2 (1238/ 675-1441/ 806)  26.6 W16-W28 S22-S13
     18  969  444 (S43W05 L= 67)  2651    493        9 ( 746/ 400-1192/ 490)  11.6 E13-W21 S45-S37
  

    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