Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  975 1773 (N38E30 L= 98)  2489    455        7 ( 853/1654-1054/1868) -22.9 E45-E19 N30-N43
      2  453 1678 (N18E61 L= 67)  3076    427        5 ( 389/1615- 529/1762)  56.3 E68-E54 N13-N24
      3  538 1666 (N19E53 L= 76)   857    121        2 ( 492/1630- 582/1711)  66.8 E55-E49 N15-N22
      4 1465 1651 (N48W10 L=139)   381     62        1 (1447/1621-1490/1675)  78.7 W08-W12 N44-N49
      5  523 1551 (N11E47 L= 81)  1014    205        4 ( 452/1473- 566/1604) -17.1 E54-E41 N07-N14
      6  596 1531 (N12E41 L= 87)   479     82        2 ( 572/1486- 623/1568) -80.5 E42-E38 N08-N14
      7  246 1500 (N03E75 L= 54)   942    163        3 ( 222/1466- 269/1532)  68.2 E80-E71 N00-N05
      8 1355 1319 (N20W11 L=139)   640    122        2 (1310/1272-1410/1351)  -7.1 W07-W15 N17-N22
      9 1519  503 (S19W48 L=177)   351     88        0 (1475/ 482-1559/ 520)  51.3 W46-W50 S21-S16
     10 1570  440 (S20W56 L=184)   573    129        1 (1521/ 404-1600/ 462)  38.7 W52-W59 S22-S18
     11 1143  297 (S44W38 L=166)  1329    255        3 (1037/ 269-1255/ 338)  23.6 W28-W45 S47-S39
  

    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