Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  653 1615 (N20E42 L=177)   284     73        0 ( 611/1602- 685/1627)   0.0 E45-E38 N19-N21
      2 1350 1598 (N40W04 L=223)   392    100        0 (1295/1580-1401/1618)  14.0 W00-W08 N38-N41
      3 1270 1411 (N24W04 L=223)   337    115        0 (1239/1365-1311/1455)  76.0 W02-W05 N20-N28
      4  699 1117 (S07E21 L=198)   198     83        1 ( 663/1091- 738/1144)   0.0 E23-E18 S08-S06
      5  221 1056 (S22E58 L=161)   722     88        2 ( 196/1006- 239/1098)  90.0 E60-E55 S24-S18
      6  168 1007 (S26E65 L=154)   448    112        2 ( 155/ 962- 186/1057)  71.6 E67-E62 S28-S22
      7  452  899 (S26E33 L=186)   314     89        0 ( 414/ 875- 494/ 927)  -8.1 E36-E29 S26-S24
      8  743  769 (S26E10 L=209)   591    212        0 ( 702/ 710- 809/ 859) -36.9 E13-E04 S28-S20
      9  317  568 (S52E51 L=168)   621     70        1 ( 291/ 543- 335/ 610)  78.7 E57-E47 S54-S48
  

    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