Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1063 1699 (N38E10 L=249)   364     86        1 (1020/1667-1099/1729) -15.9 E13-E06 N35-N39
      2 1587 1472 (N32W32 L=291)   372     88        0 (1557/1430-1617/1509) -38.7 W29-W34 N29-N33
      3 1771 1247 (N21W48 L=307)   584    152        0 (1732/1178-1790/1305) -54.5 W44-W50 N16-N23
      4  208 1162 (S07E60 L=199)   551     47        2 ( 184/1134- 229/1193) -63.4 E63-E57 S09-S05
      5  539 1077 (S10E31 L=228)   595    147        1 ( 498/1021- 599/1134)  58.0 E33-E27 S13-S05
      6  222  923 (S22E58 L=201)   867    156        0 ( 193/ 855- 244/ 995)  76.0 E61-E54 S25-S17
      7  436  882 (S22E37 L=222)   266     84        0 ( 402/ 851- 474/ 904)  38.7 E39-E34 S24-S20
      8  799  718 (S28E10 L=249)  1669    519        2 ( 616/ 622-1041/ 777)  29.1 E22-W05 S36-S20
      9 1717  654 (S15W56 L=314)   560     69        1 (1684/ 632-1747/ 681)  45.0 W52-W58 S16-S12
     10 1527  612 (S21W41 L=300)   584    152        0 (1453/ 596-1600/ 630)  21.8 W35-W46 S23-S18
     11 1679  431 (S28W68 L=327)  1973    230        2 (1599/ 381-1777/ 513)  47.5 W61-W74 S32-S20
     12  961  230 (S61W17 L=276)  1155    115        1 ( 921/ 187-1005/ 289) -26.6 W09-W26 S65-S56
  

    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