Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1409 1785 (N38W63 L=152)  1482    172        3 (1341/1728-1481/1818)  -5.0 W51-W74 N34-N39
      2 1371 1665 (N33W47 L=136)   942    163        0 (1314/1603-1421/1721)  78.7 W44-W50 N27-N37
      3 1512 1530 (N21W48 L=137)   292     62        0 (1480/1511-1545/1549) -45.0 W45-W49 N18-N22
      4 1875  960 (S22W59 L=148)   888    111        2 (1821/ 934-1910/ 994)  -4.4 W51-W64 S24-S20
      5 1601  881 (S17W31 L=120)   820    197        0 (1576/ 788-1615/ 980)  59.0 W27-W33 S22-S11
      6  236  801 (N11E61 L= 28)  1997    308        5 ( 183/ 668- 290/ 927) -76.8 E64-E57 N02-N19
      7  930  410 (S27E24 L= 65)   529     87        0 ( 884/ 374- 967/ 446)  20.6 E28-E20 S28-S24
      8 1397  331 (S49W05 L= 94)   354     87        0 (1362/ 311-1441/ 359)   0.0 W01-W09 S50-S47
      9 1278  251 (S52E09 L= 80)   398     76        0 (1241/ 229-1316/ 271)   0.0 E13-E05 S53-S50
     10  777  244 (S35E45 L= 43)   412     58        0 ( 750/ 224- 808/ 268)  14.0 E49-E41 S36-S33
  

    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