Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1354 1571 (N37W04 L= 82)  2242    532        0 (1275/1352-1467/1786) -50.6 E07-W17 N25-N54
      2  362 1633 (N14E69 L=  8)  2962    124        3 ( 311/1582- 419/1682)  14.0 E81-E61 N10-N17
      3  708 1527 (N15E33 L= 44)   385    111        0 ( 667/1498- 737/1552)  11.3 E35-E30 N12-N16
      4 1154 1500 (N25E06 L= 71)   362    111        0 (1103/1463-1195/1529)  -8.1 E09-E02 N24-N26
      5  442 1442 (N04E48 L= 29)   289     48        1 ( 413/1426- 468/1459)  45.0 E50-E45 N02-N05
      6 1599 1306 (N27W28 L=105)   273     63        1 (1579/1282-1625/1336) -26.6 W25-W30 N26-N28
      7  389 1252 (S08E45 L= 32)   321     78        0 ( 375/1206- 401/1288) -78.7 E47-E43 S10-S05
      8 1680 1204 (N24W36 L=113)   528    121        0 (1631/1176-1718/1237)  20.6 W31-W39 N21-N26
      9  934 1105 (S03E07 L= 70)   224     32        0 ( 915/1080- 953/1131) -56.3 E08-E05 S04-S01
     10  104 1039 (S25E77 L=  0)  1032     63        0 (  76/1011- 138/1059)   7.1 E89-E69 S26-S24
     11  906  959 (S12E05 L= 72)   440    170        1 ( 833/ 920- 954/1015)  63.4 E08-E02 S15-S07
     12  493  987 (S21E31 L= 46)   247     71        0 ( 456/ 960- 516/1008) -11.3 E34-E29 S21-S19
     13 1574  793 (S02W39 L=116)  1343    285        4 (1538/ 649-1612/ 915) -52.4 W33-W45 S09-N04
     14 1762  837 (N06W52 L=129)   613     94        2 (1734/ 786-1788/ 891) -71.6 W50-W54 N03-N09
     15 1806  801 (N06W58 L=135)   575     69        1 (1787/ 770-1835/ 831)  90.0 W55-W61 N03-N07
     16 1798  729 (N02W60 L=137)   537     70        2 (1765/ 705-1825/ 757)  36.9 W56-W63 N00-N03
  

    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