Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  957 1432 (N28W05 L=335)   330    104        0 ( 910/1398-1005/1465) -56.3 W03-W07 N24-N30
      2 1813 1261 (S02W58 L= 27)  3887    456        4 (1656/1080-1876/1363) -86.6 W45-W66 S13-N04
      3  913 1317 (N22W00 L=329)   199     87        0 ( 869/1307- 955/1331) -26.6 E02-W02 N20-N22
      4 1321 1291 (N11W23 L=352)   312     83        0 (1273/1269-1365/1311) -38.7 W20-W25 N09-N13
      5 1338  696 (S23W11 L=340)  2835    832        1 ( 988/ 514-1618/ 889)   9.9 E14-W32 S26-S17
      6  324  850 (N07E48 L=281)   299     44        0 ( 300/ 823- 342/ 873) -71.6 E50-E46 N05-N08
      7 1067  796 (S11E03 L=327)   425    165        0 ( 988/ 759-1136/ 840)   0.0 E07-W02 S11-S09
      8  231  786 (N05E59 L=270)  1169    186        3 ( 164/ 745- 319/ 837) -18.4 E66-E51 N01-N08
      9  145  765 (N05E72 L=257)   776     66        0 ( 122/ 731- 172/ 795)  15.3 E78-E66 N03-N06
     10  940  438 (S30E20 L=309)   359     97        0 ( 920/ 389- 972/ 478)  35.5 E23-E16 S32-S27
     11  535  306 (S29E58 L=271)  1520    188        3 ( 464/ 250- 609/ 356)   2.6 E74-E48 S32-S26
  

    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