Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  525 1290 (N16E31 L= 94)   413    120        1 ( 475/1258- 573/1313) -29.7 E34-E27 N14-N18
      2  240 1230 (N15E55 L= 70)   393    104        1 ( 212/1199- 265/1258) -38.7 E57-E51 N12-N16
      3 1230 1161 (N05W13 L=138)   180     66        0 (1197/1129-1252/1183) -53.1 W11-W14 N02-N06
      4 1694 1021 (S05W44 L=169)   368     74        0 (1662/ 978-1721/1055) -59.0 W41-W46 S07-S02
      5 1530 1018 (S05W32 L=156)   260     57        0 (1512/ 980-1547/1051) -63.4 W30-W32 S06-S02
      6  124  950 (S01E70 L= 55)   899    251        0 (  78/ 923- 170/ 980) -10.0 E80-E63 S02-N01
      7 1853  669 (S27W63 L=188)   672     59        0 (1832/ 640-1876/ 704)  63.4 W60-W65 S28-S24
      8  514  636 (S23E37 L= 88)   238     42        0 ( 497/ 610- 528/ 663) -90.0 E37-E35 S24-S21
      9  395  606 (S24E47 L= 77)   306    113        0 ( 369/ 581- 438/ 656) -68.2 E48-E44 S25-S20
     10 1094  496 (S36W02 L=127)   527    122        0 (1034/ 471-1169/ 512)   0.0 E02-W07 S37-S34
  

    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