Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1450 1736 (N56W35 L= 87)   479     50        1 (1420/1717-1481/1753)   5.7 W30-W40 N53-N58
      2  340 1599 (N29E70 L=342)   786     89        1 ( 319/1552- 354/1643) -35.0 E76-E65 N25-N32
      3  590 1504 (N26E37 L= 15)   910    260        1 ( 520/1436- 656/1539)  -9.5 E41-E29 N22-N28
      4  411 1446 (N20E50 L=  2)   394     75        3 ( 365/1423- 448/1468)  -5.7 E55-E45 N19-N21
      5 1105 1141 (N09W04 L= 56)   220     74        1 (1080/1100-1127/1176) -45.0 W01-W05 N06-N10
      6  659 1130 (N04E24 L= 28)   505    171        1 ( 586/1092- 727/1157)   5.7 E28-E18 N01-N05
      7 1711  901 (N00W48 L=100)  1037    187        0 (1653/ 821-1748/ 981)  60.9 W43-W50 S04-N05
      8  150  890 (S17E69 L=343)   674     63        0 ( 113/ 872- 180/ 910)  12.1 E77-E63 S18-S15
      9 1857  642 (S14W76 L=128)  2017    105        2 (1837/ 594-1888/ 697) -66.8 W72-W83 S17-S10
     10 1778  630 (S16W63 L=115)   429     40        0 (1751/ 610-1798/ 652)  36.9 W60-W65 S17-S14
     11 1495  284 (S43W59 L=111)   662     75        1 (1462/ 258-1511/ 315)  71.6 W54-W65 S45-S40
  

    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