Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  828 1334 (N12E18 L=104)   346    121        0 ( 794/1309- 883/1376) -15.9 E21-E14 N11-N14
      2  720 1127 (S02E20 L=102)   447    149        2 ( 676/1056- 776/1193) -29.1 E24-E15 S04-N01
      3  278 1096 (S12E52 L= 69)   573     55        1 ( 260/1055- 293/1146) -78.7 E54-E50 S14-S09
      4  965 1025 (S03E03 L=118)   304     88        2 ( 948/ 977- 981/1071) -68.2 E04-E02 S05-N00
      5  711  878 (S17E16 L=106)   989    333        2 ( 563/ 795- 813/ 990)  50.9 E24-E11 S24-S08
      6  659  722 (S27E17 L=105)   346    115        0 ( 625/ 676- 704/ 764)  63.4 E19-E14 S31-S23
      7 1344  564 (S22W31 L=153)   270     70        0 (1302/ 546-1379/ 583)  36.9 W28-W33 S23-S20
      8 1724  525 (S16W64 L=186)   478     72        0 (1694/ 499-1752/ 549)  45.0 W61-W66 S17-S13
      9 1241  357 (S38W34 L=156)  1104    273        1 (1112/ 335-1374/ 377)  21.3 W24-W43 S42-S33
  

    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