Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  932 1619 (N37E17 L=354)   439    145        0 ( 887/1570- 988/1673) -28.6 E22-E11 N34-N40
      2 1269 1548 (N37W10 L= 21)   299     54        1 (1248/1531-1304/1569) -14.0 W08-W12 N35-N39
      3  507 1378 (N15E40 L=331)   646    197        1 ( 477/1277- 529/1462) -47.7 E45-E35 N09-N20
      4 1736 1444 (N37W55 L= 66)   469     56        1 (1697/1428-1763/1459)  11.3 W49-W59 N36-N38
      5  237 1371 (N11E67 L=304)   921     88        2 ( 218/1332- 256/1424) -63.4 E70-E63 N08-N14
      6  752 1225 (N09E19 L=352)   867    221        0 ( 709/1118- 795/1315)  80.5 E20-E17 N02-N14
      7  975 1177 (N09E05 L=  6)   312    137        1 ( 940/1122-1011/1232) -45.0 E07-E02 N06-N11
      8  435 1117 (S02E40 L=331)   395     84        0 ( 407/1088- 472/1156) -31.0 E42-E36 S03-N00
      9  968  663 (S23W01 L= 12)  2122    658        3 ( 718/ 565-1239/ 774)  34.1 E15-W16 S32-S11
     10  333  601 (S37E52 L=319)  1014    204        0 ( 287/ 526- 396/ 695) -55.0 E55-E46 S41-S30
     11  481  630 (S32E34 L=337)   404     90        0 ( 441/ 619- 531/ 648)   7.1 E38-E30 S32-S31
  

    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