Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  729 1740 (N61W05 L= 19)   550    112        1 ( 700/1716- 770/1768) -56.3 W01-W08 N57-N63
      2 1127 1710 (N43W34 L= 47)   803    154        2 (1051/1667-1194/1742) -60.9 W29-W36 N37-N46
      3 1494 1688 (N28W58 L= 72)   720    148        0 (1413/1670-1559/1702) -32.5 W51-W63 N24-N31
      4  573 1633 (N58E17 L=357)   639     98        0 ( 536/1595- 612/1663)  33.7 E21-E12 N55-N60
      5  388 1547 (N57E42 L=331)   825    173        0 ( 333/1491- 459/1593)  16.4 E49-E32 N53-N58
      6 1183 1523 (N29W26 L= 39)   299     12        0 (1160/1504-1206/1542) -18.4 W24-W27 N27-N30
      7 1465 1464 (N17W40 L= 54)   329     78        3 (1433/1430-1494/1491) -78.7 W37-W42 N14-N19
      8 1252 1342 (N17W22 L= 36)   413    212        0 (1150/1325-1355/1368) -38.7 W16-W27 N13-N21
      9  308 1123 (N29E47 L=327)   364     86        0 ( 266/1105- 338/1154)   8.1 E51-E43 N28-N30
     10 1476  627 (S30W17 L= 30)   672    145        0 (1409/ 591-1529/ 674)   5.7 W11-W21 S31-S28
     11 1252  326 (S44E08 L=  6)   992    187        2 (1171/ 275-1325/ 380)   7.1 E16-W00 S45-S41
     12  906  225 (S40E39 L=335)   665    124        0 ( 876/ 184- 940/ 284)  63.4 E42-E34 S44-S34
  

    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