Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  836 1797 (N49E22 L=344)   393     54        1 ( 809/1771- 861/1823)  80.5 E24-E19 N45-N51
      2  233  732 (S23E61 L=305)  1340    178        7 ( 195/ 678- 264/ 807)  76.0 E67-E57 S26-S18
      3  289  759 (S21E53 L=313)   403     68        1 ( 267/ 728- 312/ 797) -90.0 E55-E51 S23-S19
      4 1705  649 (S23W53 L= 59)   306     36        0 (1682/ 632-1723/ 670)  45.0 W50-W54 S23-S21
      5 1469  581 (S29W34 L= 40)   632    160        1 (1414/ 540-1525/ 626)  45.0 W30-W38 S32-S25
      6 1325  531 (S33W24 L= 31)   284     25        0 (1301/ 509-1346/ 549) -14.0 W22-W26 S34-S32
      7  292  515 (S37E66 L=300)   616     43        1 ( 268/ 492- 321/ 541)  14.0 E72-E60 S38-S35
      8 1184  430 (S41W16 L= 22)   227     72        0 (1148/ 419-1224/ 444)   9.5 W12-W19 S42-S40
  

    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