Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  552 1635 (N17E50 L=170)  1365    232        1 ( 528/1523- 599/1725) -45.0 E57-E44 N09-N22
      2 1535 1667 (N52W17 L=236)   795    146        0 (1476/1632-1591/1704)  -4.4 W10-W23 N49-N53
      3  174 1180 (S17E66 L=154)  1362    160        2 ( 148/1098- 209/1246)  68.2 E68-E61 S22-S12
      4  383 1061 (S20E41 L=178)   955    187        2 ( 303/1021- 468/1090)  24.8 E47-E34 S23-S16
      5 1720  902 (N08W46 L=265)   300     52        0 (1695/ 879-1744/ 922)   0.0 W43-W48 N07-N08
      6 1526  579 (S15W45 L=264)   286     80        1 (1492/ 557-1565/ 593)  36.9 W41-W47 S16-S13
      7 1627  525 (S14W55 L=275)  1270    246        1 (1565/ 455-1679/ 578)  90.0 W50-W60 S19-S10
      8 1523  473 (S20W51 L=270)   439     84        1 (1494/ 435-1548/ 511) -71.6 W48-W52 S23-S17
      9  858  396 (S48W12 L=231)  1286    256        0 ( 747/ 339- 974/ 441)   7.8 W01-W23 S50-S44
     10 1194  191 (S48W56 L=275)  1362    218        2 (1097/ 133-1268/ 235)   1.9 W42-W74 S50-S45
  

    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