Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1549 1640 (N29W53 L=270)   342     48        0 (1523/1617-1569/1660) -56.3 W50-W54 N27-N30
      2 1142 1639 (N37W21 L=238)   264     33        0 (1113/1625-1166/1655)   0.0 W18-W22 N35-N37
      3  634 1362 (N27E20 L=198)  2140    536        2 ( 461/1245- 852/1456)  13.6 E33-E03 N22-N29
      4  296 1328 (N31E49 L=168)  1836    217        2 ( 192/1297- 405/1370) -14.5 E64-E37 N27-N34
      5 1746 1310 (N06W54 L=271)   337     41        0 (1716/1297-1773/1323)   0.0 W50-W56 N04-N06
      6 1815 1044 (S11W55 L=272)   626    103        2 (1787/ 998-1849/1091) -63.4 W52-W58 S13-S07
      7 1743  927 (S17W47 L=264)   312     16        1 (1721/ 906-1764/ 944) -14.0 W44-W48 S17-S15
      8 1326  836 (S15W15 L=233)   242     59        0 (1296/ 815-1358/ 856)  14.0 W13-W17 S15-S14
      9 1507  750 (S23W27 L=244)   437    101        0 (1492/ 697-1524/ 806) -81.9 W25-W27 S26-S19
     10  766  765 (S10E20 L=198)   196     54        0 ( 753/ 732- 781/ 796)  63.4 E20-E18 S12-S08
     11 1607  717 (S27W34 L=252)   539    102        0 (1560/ 688-1658/ 754)  14.0 W30-W38 S28-S25
     12  444  499 (S21E51 L=166)   432    109        0 ( 393/ 477- 481/ 537)  11.3 E56-E46 S22-S19
  

    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