Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  489 1231 (N31E29 L= 37)  1845    351        1 ( 349/1180- 638/1269) -14.6 E41-E18 N25-N34
      2  415 1013 (N20E38 L= 28)  1185    321        2 ( 291/ 955- 550/1090) -35.8 E48-E29 N13-N26
      3  684  930 (N09E22 L= 44)  1048    268        0 ( 636/ 869- 757/1034)  78.7 E24-E17 N04-N14
      4  183  943 (N20E63 L=  3)  1223    193        3 ( 139/ 863- 242/1022) -47.7 E68-E57 N14-N25
      5 1818  893 (S26W50 L=116)   394     54        1 (1794/ 870-1839/ 923)  11.3 W47-W52 S26-S24
      6  409  896 (N13E41 L= 25)   293     59        3 ( 383/ 873- 451/ 919) -36.9 E42-E38 N11-N14
      7  849  769 (S04E16 L= 49)   188     56        0 ( 821/ 749- 869/ 796)  14.0 E18-E14 S05-S02
      8 1684  566 (S43W35 L=101)  1291    247        3 (1579/ 506-1782/ 626)   0.0 W23-W47 S44-S40
      9  433  383 (S18E63 L=  3)  3347    445        7 ( 320/ 255- 542/ 499)  84.3 E71-E57 S28-S08
     10  918  280 (S37E33 L= 33)  1907    221        5 ( 829/ 198- 993/ 357)  20.0 E45-E23 S41-S32
  

    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