Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1425 1466 (N24W37 L=199)  1519    284        2 (1308/1372-1502/1556) -68.2 W32-W40 N17-N32
      2 1624 1510 (N22W56 L=217)   471    141        1 (1559/1483-1676/1527) -29.1 W49-W60 N18-N23
      3 1020 1222 (N19W04 L=165)   193     94        0 ( 980/1188-1054/1260) -68.2 W01-W05 N16-N21
      4 1282 1188 (N11W19 L=180)   514    278        0 (1226/1125-1372/1225)  11.3 W14-W24 N07-N13
      5 1639  929 (S12W37 L=199)   332    160        0 (1607/ 895-1670/ 965)  31.0 W34-W40 S13-S09
      6 1400  942 (S05W21 L=182)   357    101        2 (1355/ 925-1458/ 961) -33.7 W18-W24 S07-S03
      7 1664  748 (S24W38 L=199)  1748    492        0 (1519/ 596-1802/ 931)  27.4 W25-W52 S30-S15
      8 1511  850 (S13W27 L=188)   379    117        0 (1482/ 806-1528/ 901)  68.2 W24-W28 S15-S10
      9  453  870 (N07E40 L=122)   327     85        0 ( 415/ 851- 492/ 886)   0.0 E42-E36 N05-N07
     10 1745  691 (S30W47 L=208)  1198    243        1 (1675/ 597-1812/ 774)  32.7 W39-W54 S35-S26
     11  638  576 (S14E34 L=128)   243     66        2 ( 604/ 559- 668/ 589)   0.0 E36-E31 S15-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