Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1689 1524 (N20W61 L= 50)   637    120        2 (1654/1475-1730/1576) -74.1 W58-W63 N16-N23
      2 1783 1389 (N10W64 L= 53)   491     64        1 (1767/1363-1803/1415) -90.0 W61-W66 N08-N11
      3 1650 1362 (N12W49 L= 39)   551    116        2 (1606/1324-1694/1391) -51.3 W45-W52 N09-N14
      4  672 1148 (N20E20 L=330)  1010    313        4 ( 534/1082- 784/1228) -40.9 E28-E12 N14-N27
      5  276 1144 (N26E51 L=298)   799    154        3 ( 243/1065- 300/1192) -66.0 E54-E48 N20-N29
      6 1015 1151 (N15W02 L=351)   402    156        2 ( 946/1100-1100/1181) -40.6 E02-W06 N10-N17
      7  220  824 (N07E62 L=287)  1323    303        1 ( 162/ 695- 283/ 952) -68.2 E67-E58 S01-N14
      8  809  765 (S05E18 L=331)  1514    632        7 ( 693/ 638- 947/ 923)  34.5 E25-E06 S12-N03
      9  423  456 (S19E56 L=293)  1518    387        1 ( 408/ 347- 444/ 567)  53.1 E64-E50 S27-S11
  

    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