Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1146 1638 (N41W27 L=338)   282     52        1 (1125/1608-1162/1660)  76.0 W24-W28 N38-N42
      2 1737 1465 (N14W63 L= 14)   877    121        1 (1709/1430-1773/1532)  81.9 W59-W65 N10-N17
      3 1201 1495 (N30W24 L=335)   290     62        1 (1183/1467-1235/1519) -56.3 W21-W25 N29-N32
      4 1786 1410 (N09W65 L= 16)   723     87        1 (1763/1381-1812/1435)   8.1 W62-W70 N07-N10
      5 1509 1220 (N07W34 L=345)   806    225        1 (1488/1103-1532/1326)  69.0 W30-W37 N00-N13
      6  155  973 (N18E66 L=246)  1683    362        4 (  93/ 898- 219/1047) -39.3 E72-E59 N13-N22
      7 1617  790 (S21W32 L=344)   671    186        1 (1555/ 710-1661/ 854)  35.0 W26-W36 S24-S17
      8  317  721 (N01E55 L=256)   517     77        0 ( 291/ 679- 353/ 759) -59.0 E56-E52 S01-N04
      9 1379  452 (S37W11 L=322)   342     36        1 (1355/ 432-1403/ 469)   0.0 W08-W12 S39-S35
  

    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