Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1180 1758 (N44E17 L=197)   644    161        0 (1113/1696-1228/1820) -21.8 E25-E10 N40-N47
      2  344 1406 (S02E56 L=157)  1441    221        3 ( 306/1298- 392/1482) -34.7 E63-E48 S07-N03
      3 1192 1371 (N18E00 L=213)   796    192        0 (1115/1290-1261/1424) -15.3 E06-W05 N14-N19
      4  218 1065 (S24E59 L=154)   405     45        2 ( 201/1036- 236/1096) -56.3 E60-E56 S25-S21
      5  783  905 (S19E10 L=203)  1101    437        0 ( 595/ 813- 966/ 962)  41.8 E20-E01 S29-S11
      6  403  798 (S36E36 L=177)   734    174        1 ( 371/ 725- 464/ 864)  65.6 E39-E31 S41-S30
      7  912  748 (S25W02 L=215)   604    178        1 ( 840/ 707-1001/ 776)  48.8 E01-W06 S29-S20
      8  978  621 (S30W10 L=223)  1007    288        1 ( 844/ 602-1110/ 647)  31.0 W01-W17 S34-S24
      9 1266  496 (S28W33 L=246)   330     88        0 (1224/ 476-1307/ 521)   9.5 W29-W35 S28-S26
     10  418  445 (S62E39 L=174)   505     37        0 ( 397/ 426- 442/ 472) -26.6 E44-E34 S63-S59
     11 1542  439 (S21W56 L=269)   639     88        1 (1497/ 413-1587/ 467)  63.4 W53-W58 S24-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