Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  958 1901 (N62E44 L=111)   699     46        1 ( 939/1876- 978/1928) -14.7 E56-E36 N59-N65
      2  995 1867 (N59E33 L=123)   485     46        1 ( 972/1848-1019/1893)  63.4 E38-E29 N56-N62
      3  652 1578 (N27E38 L=117)   482    165        1 ( 609/1553- 716/1627)  -5.7 E43-E33 N25-N29
      4  323 1549 (N20E69 L= 86)  1206     87        2 ( 278/1522- 368/1569)   0.0 E78-E61 N18-N20
      5  554 1465 (N18E40 L=115)   415     80        1 ( 534/1436- 579/1492) -31.0 E42-E37 N16-N19
      6 1846 1314 (N32W65 L=220)  1792    261        2 (1759/1242-1910/1397) -16.9 W53-W79 N27-N37
      7 1565 1266 (N23W32 L=187)   305    103        0 (1554/1215-1579/1313) -80.5 W31-W33 N20-N26
      8 1819  919 (N07W59 L=214)   462    117        0 (1794/ 868-1843/ 966)  60.3 W57-W61 N03-N10
      9  236  850 (S25E57 L= 98)   762    127        0 ( 214/ 795- 272/ 904)  48.8 E60-E52 S28-S20
     10  375  801 (S25E42 L=113)   319     94        0 ( 340/ 786- 418/ 831)  -7.1 E45-E37 S26-S23
     11 1038  606 (S26W09 L=164)  2111    549        7 ( 818/ 509-1250/ 732)  42.4 E04-W19 S35-S14
     12 1543  618 (S16W43 L=198)   656    172        1 (1489/ 574-1611/ 656)  45.0 W39-W47 S19-S12
  

    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