Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1383 1556 (N32W36 L= 32)   474    168        0 (1313/1538-1451/1581)  -4.8 W30-W42 N30-N33
      2 1082 1428 (N30W11 L=  7)   728    175        2 (1009/1375-1173/1465)   9.5 W05-W17 N27-N31
      3 1648 1348 (N14W48 L= 44)  1201    183        1 (1575/1283-1718/1384) -45.0 W42-W52 N08-N17
      4 1110  949 (N01W04 L=360)  1846    456        7 (1006/ 757-1159/1165) -87.8 W00-W06 S11-N15
      5  506 1063 (N16E33 L=323)   456    118        0 ( 457/1043- 558/1093) -23.2 E36-E29 N14-N17
      6  282  979 (N12E52 L=304)   309     50        1 ( 263/ 948- 296/1004) -90.0 E54-E51 N10-N13
      7 1415  895 (S08W22 L= 18)   779    215        2 (1334/ 834-1508/ 944) -45.0 W17-W27 S13-S03
      8 1868  865 (S20W60 L= 56)   466     81        2 (1850/ 830-1881/ 903)  63.4 W57-W61 S21-S17
      9 1284  788 (S12W12 L=  8)   391    106        1 (1265/ 732-1309/ 844)  66.8 W10-W14 S15-S08
     10 1295  583 (S26W10 L=  6)   627    137        0 (1269/ 518-1328/ 650)  58.0 W08-W13 S29-S21
     11  951  345 (S37E18 L=338)  1184    307        1 ( 876/ 276- 998/ 461)  42.9 E26-E12 S42-S28
  

    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