Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1110 1877 (N58E31 L=219)   991     87        0 (1086/1833-1134/1918) -42.7 E40-E26 N52-N64
      2  753 1733 (N33E42 L=208)   388     77        0 ( 736/1701- 770/1776) -40.6 E46-E39 N30-N36
      3 1276 1662 (N44E02 L=248)   619     81        2 (1244/1622-1305/1702) -45.0 E05-W01 N40-N46
      4 1620 1554 (N48W32 L=282)  1278    217        3 (1542/1472-1698/1609) -15.5 W23-W41 N42-N50
      5 1626 1254 (N26W33 L=283)  2080    512        5 (1478/1053-1768/1426) -27.3 W19-W48 N17-N32
      6  178 1004 (S24E64 L=186)  1785    193        1 ( 143/ 908- 214/1081)  65.6 E70-E58 S29-S18
      7 1099  612 (S25W16 L=266)  1261    292        3 ( 972/ 563-1236/ 652)  40.2 W09-W23 S31-S19
      8 1316  527 (S24W33 L=283)   361    154        0 (1264/ 483-1362/ 586)  77.5 W31-W33 S28-S19
      9 1426  346 (S32W53 L=303)  2099    314        4 (1260/ 321-1544/ 381)  23.5 W40-W64 S37-S27
     10  976  160 (S63W44 L=294)  1692    154        2 ( 943/ 102-1027/ 224) -13.1 W30-W73 S66-S56
  

    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