Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1196 1704 (N40W29 L=110)   623    102        0 (1158/1663-1235/1738)  18.4 W24-W33 N37-N41
      2  814 1653 (N45E01 L= 80)   272     39        1 ( 792/1634- 840/1675) -71.6 E02-W00 N43-N46
      3 1634 1581 (N22W55 L=137)   789    145        1 (1590/1543-1671/1606)   6.3 W50-W59 N19-N23
      4  913 1566 (N36W04 L= 85)   266      7        0 ( 890/1550- 936/1582) -33.7 W02-W05 N35-N37
      5  794 1344 (N24E08 L= 73)   241     61        0 ( 754/1328- 831/1361) -21.8 E10-E05 N23-N25
      6 1675 1035 (S10W41 L=123)  1947    427        0 (1576/ 864-1823/1164)  25.6 W31-W56 S18-S04
      7 1874 1079 (S11W61 L=142)   653    102        0 (1839/1037-1904/1112) -45.0 W57-W63 S14-S09
      8  795  289 (S40E36 L= 46)   468     97        0 ( 778/ 256- 817/ 326)  50.2 E38-E32 S43-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