Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  933 1741 (N51W21 L=142)   522     87        0 ( 907/1700- 969/1775)  24.0 W16-W25 N48-N53
      2 1300 1717 (N36W44 L=166)   534    115        1 (1244/1682-1359/1743) -11.3 W38-W49 N34-N37
      3 1261 1603 (N31W33 L=155)   305     54        0 (1229/1589-1294/1617) -21.8 W30-W35 N29-N31
      4  818 1361 (N29E03 L=118)   221     29        0 ( 793/1341- 839/1381) -45.0 E04-E01 N27-N30
      5  249 1348 (N43E54 L= 67)   435     60        0 ( 226/1314- 275/1377)  36.9 E57-E50 N40-N44
      6 1278 1333 (N15W23 L=144)   261     92        0 (1232/1306-1315/1350) -45.0 W20-W24 N12-N16
      7 1743 1278 (S01W52 L=174)   954    140        2 (1695/1241-1805/1313) -39.8 W47-W57 S04-N01
      8  338 1096 (N25E43 L= 79)  2532    504        3 ( 194/ 999- 528/1283) -38.7 E60-E29 N15-N39
      9 1819  949 (S21W50 L=172)   500    110        0 (1794/ 898-1834/1007)  51.3 W47-W52 S23-S18
     10 1922  944 (S25W65 L=186)   734     60        0 (1891/ 926-1955/ 966)   0.0 W59-W72 S26-S23
     11  255  809 (N10E56 L= 65)  1051    138        3 ( 214/ 746- 303/ 846)  50.2 E61-E51 N06-N12
     12 1288  745 (S17W08 L=129)   687    312        0 (1187/ 653-1403/ 840)  17.4 W00-W16 S20-S14
     13  911  773 (S06E12 L=109)  1170    342        0 ( 843/ 659-1040/ 815)   0.0 E18-E05 S11-S02
     14 1602  655 (S32W27 L=149)  1421    340        0 (1511/ 562-1715/ 795)  24.2 W18-W38 S36-S26
     15  288  713 (N04E57 L= 65)  1643    215        4 ( 240/ 631- 359/ 783)  26.6 E65-E49 N00-N07
     16  201  675 (N03E69 L= 53)  1160    123        4 ( 157/ 617- 243/ 717) -56.3 E74-E62 N00-N06
     17  481  594 (S06E45 L= 77)   346     61        0 ( 452/ 577- 515/ 613) -11.3 E47-E41 S07-S05
     18  840  163 (S44E50 L= 72)   596     92        0 ( 799/ 143- 875/ 196)   4.8 E55-E42 S45-S42
  

    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