Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1452 1536 (N38W32 L=342)   284     40        1 (1439/1507-1473/1562)  76.0 W29-W33 N35-N39
      2  722 1363 (N22E22 L=288)   286    101        0 ( 679/1341- 765/1391) -26.6 E25-E19 N20-N23
      3 1359 1247 (N18W21 L=331)   994    438        0 (1193/1106-1503/1348) -33.7 W09-W30 N09-N23
      4 1575 1246 (N19W37 L=347)   251     80        0 (1552/1216-1608/1265) -14.0 W35-W39 N16-N19
      5 1667  800 (S09W46 L=356)   361     65        2 (1628/ 777-1700/ 824)  -9.5 W42-W48 S10-S07
      6  876  694 (S19E08 L=302)   510    157        2 ( 835/ 636- 920/ 761)  58.0 E10-E05 S23-S14
      7 1755  716 (S14W57 L=  7)   365     56        0 (1739/ 689-1772/ 743)  71.6 W54-W58 S15-S12
      8  446  506 (S36E45 L=265)   577     82        0 ( 410/ 488- 491/ 529)  -7.1 E48-E40 S36-S33
      9 1206  393 (S39W19 L=329)   324     42        1 (1179/ 370-1227/ 417)  18.4 W16-W20 S40-S36
  

    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