Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  699 1337 (N30E17 L=  4)  1436    336        5 ( 551/1280- 846/1395)  -2.6 E28-E06 N26-N32
      2 1327 1228 (N14W22 L= 44)   588    132        2 (1262/1195-1388/1259) -35.5 W18-W25 N10-N15
      3  187 1045 (N16E62 L=319)  1027    239        1 ( 138/1009- 242/1094) -18.4 E68-E56 N13-N19
      4 1681  976 (S08W42 L= 64)   933    241        1 (1635/ 864-1723/1084)  47.7 W37-W47 S13-S02
      5  697  740 (S07E24 L=358)  1563    428        9 ( 612/ 588- 743/ 890)  57.1 E33-E20 S15-N02
      6  247  806 (N01E60 L=322)   393     58        2 ( 234/ 771- 264/ 844) -90.0 E61-E57 N00-N03
      7  963  806 (S06E07 L= 15)   440    138        3 ( 888/ 778-1023/ 830) -20.6 E11-E03 S08-S04
      8 1763  746 (S24W49 L= 70)  1148    126        1 (1714/ 698-1803/ 787)  29.7 W43-W52 S26-S21
      9 1584  710 (S23W33 L= 55)   376     75        2 (1567/ 672-1598/ 748) -80.5 W31-W33 S25-S19
     10  899  511 (S24E16 L=  6)   274     60        0 ( 876/ 483- 932/ 539)  21.8 E18-E13 S25-S22
     11  671  290 (S38E42 L=340)  1085    186        2 ( 630/ 253- 700/ 336)  41.2 E47-E38 S41-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