Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1120 1715 (N39E18 L= 97)  2275    473        3 ( 954/1596-1216/1820)  -7.1 E32-E07 N30-N43
      2  523 1633 (N17E52 L= 62)  2327    326        9 ( 424/1535- 627/1724)  66.8 E55-E46 N08-N24
      3  578 1601 (N16E46 L= 68)   653    129        2 ( 544/1553- 606/1650) -66.8 E48-E43 N12-N19
      4  345 1515 (N06E61 L= 53)  1512    193        3 ( 298/1442- 392/1599) -65.6 E66-E57 N00-N11
      5  706 1455 (N10E31 L= 83)  1077    337        2 ( 619/1325- 768/1542) -26.6 E38-E22 N04-N15
      6  777 1466 (N13E27 L= 87)   346    147        0 ( 749/1407- 840/1514)  82.9 E28-E24 N08-N16
      7 1551 1227 (N21W25 L=140)   871    168        2 (1500/1180-1594/1287) -26.6 W21-W29 N17-N24
      8 1223  234 (S45W48 L=162)  1639    261        2 (1103/ 208-1337/ 269)  26.6 W38-W55 S48-S39
  

    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