Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1214 1649 (N35W32 L=134)   330     74        0 (1186/1627-1263/1665)  -9.5 W29-W35 N33-N35
      2 1619 1364 (N07W46 L=148)   282     69        1 (1598/1332-1645/1385)  11.3 W43-W48 N05-N08
      3 1803 1212 (S06W56 L=158)  1315    331        4 (1721/1140-1895/1318) -49.8 W52-W64 S12-N01
      4 1584 1079 (S07W34 L=136)   872    233        1 (1553/ 949-1604/1193)  74.1 W31-W36 S14-N00
      5  828 1079 (N12E10 L= 92)   286     91        0 ( 778/1064- 877/1096)  -9.5 E12-E06 N10-N12
      6  285 1010 (N20E49 L= 53)   459    123        1 ( 241/ 994- 331/1027) -20.6 E52-E44 N18-N21
      7  424  974 (N15E38 L= 64)   400     98        1 ( 386/ 923- 452/1010) -63.4 E40-E36 N11-N17
      8 1251  821 (S13W08 L=109)   329     75        0 (1209/ 797-1294/ 839)   0.0 W04-W10 S13-S11
      9  964  572 (S20E15 L= 86)   498    226        0 ( 944/ 467- 989/ 667)  56.3 E19-E11 S26-S14
     10 1422  593 (S31W13 L=115)   475     70        2 (1394/ 561-1456/ 628)  76.0 W11-W16 S33-S29
     11  382  575 (S07E54 L= 48)  1098    187        0 ( 311/ 545- 462/ 620) -31.0 E59-E47 S09-S02
  

    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