Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1396 1636 (N46W10 L=127)   995    225        4 (1370/1570-1422/1722) -71.6 W06-W13 N40-N53
      2  497 1657 (N23E55 L= 63)   881    245        3 ( 470/1597- 528/1719) -52.1 E59-E51 N18-N27
      3  915 1661 (N33E25 L= 92)   274     25        0 ( 895/1642- 930/1684)  90.0 E26-E23 N31-N34
      4  747 1532 (N20E30 L= 87)   413    142        1 ( 717/1490- 771/1578) -71.6 E33-E27 N17-N23
      5  662 1498 (N16E34 L= 83)   253    122        0 ( 637/1465- 700/1542) -26.6 E37-E31 N14-N18
      6  988 1365 (N16E10 L=107)   362    122        3 ( 971/1294-1004/1426) -63.4 E12-E07 N11-N19
      7  237 1200 (S10E59 L= 59)   718    139        2 ( 207/1148- 270/1249) -81.9 E61-E55 S13-S06
      8 1313 1208 (N14W13 L=130)   288     60        2 (1277/1180-1343/1238)  36.9 W10-W14 N12-N16
      9  652  752 (S27E17 L=101)   227     97        0 ( 606/ 737- 687/ 778)  18.4 E20-E14 S28-S25
     10 1819  703 (S01W66 L=183)   533    117        0 (1788/ 646-1852/ 758)  74.1 W64-W67 S05-N02
     11 1344  263 (S40W50 L=168)   791    178        3 (1276/ 248-1407/ 280)  15.3 W44-W55 S42-S37
  

    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