Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  341 1507 (N55E50 L= 13)   628     70        1 ( 320/1472- 371/1534) -29.1 E56-E44 N52-N57
      2  562 1502 (N48E20 L= 43)   327     51        1 ( 533/1473- 587/1521) -26.6 E23-E17 N46-N50
      3 1876 1046 (S18W60 L=123)  1870    189        3 (1835/ 964-1910/1138) -84.8 W57-W65 S22-S11
      4  227  880 (N16E60 L=  3)  4181    517        6 ( 160/ 717- 367/1120) -59.3 E67-E49 N04-N31
      5 1736  731 (S32W40 L=104)   646     91        0 (1707/ 686-1760/ 785)  51.3 W37-W43 S34-S28
      6 1207  518 (S28E03 L= 60)   237     55        0 (1176/ 502-1233/ 540)  14.0 E05-E01 S28-S26
      7  391  437 (S13E64 L=360)   603     80        1 ( 363/ 402- 421/ 472) -78.7 E66-E61 S15-S10
      8 1046  342 (S35E20 L= 43)  1005    202        0 ( 977/ 280-1101/ 442)  31.0 E27-E12 S38-S29
      9  866  185 (S43E46 L= 17)  1124    141        1 ( 829/ 123- 907/ 240)  19.8 E61-E36 S48-S38
  

    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