Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1243 1338 (N21W20 L=195)   525    207        0 (1139/1319-1334/1362) -28.6 W13-W25 N18-N24
      2 1853  892 (S20W59 L=234)   483    107        1 (1827/ 846-1878/ 941)  33.7 W55-W62 S22-S17
      3 1720  904 (S15W44 L=219)   347     77        0 (1684/ 868-1748/ 940)  26.6 W40-W47 S16-S13
      4 1106  747 (S11E00 L=175)   210    105        0 (1058/ 717-1145/ 767) -53.1 E02-W01 S13-S09
      5 1793  669 (S33W54 L=229)   496     44        1 (1776/ 644-1810/ 704) -59.0 W51-W57 S35-S30
      6 1112  645 (S18E02 L=173)   417    179        0 (1040/ 616-1180/ 685)   0.0 E05-W03 S18-S15
      7 1566  566 (S34W30 L=204)  1599    273        0 (1502/ 461-1627/ 676)  50.2 W24-W34 S40-S27
      8  270  601 (S09E66 L=109)   726     82        2 ( 252/ 559- 285/ 644)  71.6 E68-E63 S11-S05
      9  625  524 (S17E36 L=139)  1077    296        0 ( 476/ 501- 746/ 548)  -5.7 E47-E27 S19-S14
     10 1036  373 (S35E13 L=161)  3034    524        6 ( 865/ 218-1202/ 546)  29.9 E32-W01 S46-S25
  

    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