Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  625 1539 (N12E39 L=135)   440    171        2 ( 606/1455- 643/1620) -59.0 E43-E36 N07-N17
      2 1678 1495 (N44W33 L=208)   680    119        0 (1629/1460-1728/1525)   0.0 W27-W39 N42-N45
      3  297 1402 (S03E60 L=115)   541     77        1 ( 276/1367- 317/1440) -33.7 E63-E56 S04-N00
      4  260 1172 (S17E56 L=119)   391     92        1 ( 233/1127- 282/1212) -29.7 E59-E52 S18-S14
      5  203 1048 (S25E60 L=115)   499    123        0 ( 156/1014- 248/1092)  28.6 E66-E55 S27-S21
      6  942  888 (S16E01 L=174)   289    116        0 ( 910/ 841- 982/ 935)  90.0 E02-W00 S18-S12
      7  380  846 (S33E39 L=136)   667    181        2 ( 345/ 763- 425/ 916)  63.4 E41-E35 S39-S27
      8  272  843 (S36E50 L=124)   433    127        1 ( 260/ 778- 285/ 900) -81.9 E52-E49 S39-S32
      9  489  831 (S31E29 L=146)   280     83        2 ( 465/ 796- 526/ 867) -31.0 E31-E25 S32-S29
     10  586  778 (S32E20 L=154)   540    140        1 ( 524/ 744- 639/ 824)  41.2 E25-E17 S34-S27
     11 1162  660 (S22W19 L=194)   673    281        1 (1121/ 577-1237/ 737) -79.7 W16-W24 S26-S15
     12 1377  559 (S21W36 L=211)   254     71        0 (1345/ 540-1407/ 581)  63.4 W34-W37 S22-S18
  

    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