Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  810 1533 (N24E26 L=276)   363    128        0 ( 770/1483- 839/1570) -29.7 E29-E22 N21-N26
      2 1176 1457 (N28E00 L=302)   953    296        2 (1057/1363-1288/1509) -10.6 E08-W08 N23-N29
      3 1361 1428 (N30W13 L=315)   333     83        0 (1325/1392-1396/1470) -33.7 W09-W15 N27-N32
      4 1861 1189 (N27W62 L=  4)  2362    224        3 (1819/1083-1893/1299)  70.3 W55-W67 N20-N34
      5  419 1072 (S11E40 L=262)   636    100        2 ( 374/1042- 455/1098)  33.7 E43-E36 S13-S09
      6 1738 1000 (N12W48 L=350)   663     84        0 (1716/ 953-1761/1042)  80.5 W46-W50 N09-N15
      7 1312  721 (S13W24 L=326)   279    109        0 (1267/ 710-1368/ 742)   8.1 W20-W27 S14-S12
      8 1426  573 (S20W37 L=339)  1938    361        2 (1344/ 446-1561/ 674)  70.6 W33-W43 S28-S10
      9  450  559 (S44E35 L=267)   357     47        2 ( 435/ 530- 468/ 586) -76.0 E36-E32 S46-S41
     10  897  542 (S34W03 L=305)   635    168        1 ( 813/ 513- 969/ 571)  29.1 E02-W07 S37-S30
  

    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