Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  598 1647 (N19E46 L=277)   680    114        3 ( 566/1601- 635/1701) -74.1 E47-E43 N15-N22
      2 1311 1166 (N08W12 L=334)   515    115        0 (1264/1102-1341/1222) -35.5 W08-W15 N05-N10
      3  543 1151 (S12E31 L=291)   360     88        0 ( 505/1129- 576/1169)  36.9 E33-E28 S13-S10
      4  426 1096 (S17E39 L=284)   287     76        0 ( 391/1078- 460/1112)   0.0 E41-E35 S17-S16
      5 1665 1005 (N10W39 L=  1)   367     69        0 (1650/ 963-1680/1051) -53.1 W36-W40 N08-N12
      6 1525  920 (N01W31 L=354)   752    255        0 (1452/ 806-1561/1027) -33.7 W24-W36 S04-N05
      7  821  752 (S27E05 L=318)   864    316        1 ( 671/ 712- 986/ 785)  11.3 E14-W06 S31-S24
      8 1674  674 (S07W51 L= 13)   449    171        1 (1612/ 648-1727/ 691)  21.8 W46-W56 S09-S05
      9 1599  330 (S26W66 L= 28)   578    162        1 (1574/ 313-1630/ 355)  15.9 W62-W69 S26-S24
  

    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