Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1654 1415 (N09W51 L=179)   517    156        0 (1593/1355-1715/1486) -70.0 W49-W53 N03-N14
      2 1280 1178 (N06W18 L=146)   204     70        0 (1253/1143-1310/1216)  36.9 W16-W20 N04-N07
      3 1415 1151 (N01W25 L=153)   274     69        1 (1379/1113-1452/1182) -63.4 W23-W26 S01-N03
      4  527 1119 (N21E28 L=100)   225     80        0 ( 506/1078- 557/1156)  31.0 E29-E24 N19-N22
      5 1558 1015 (S10W31 L=159)   361    141        0 (1507/1001-1625/1030) -32.0 W27-W35 S12-S07
      6  411  911 (N12E40 L= 88)   366    170        1 ( 343/ 892- 470/ 936) -11.3 E45-E35 N10-N13
      7  164  881 (N15E64 L= 64)   866    141        2 ( 131/ 810- 199/ 928) -69.4 E67-E61 N10-N18
      8 1776  847 (S26W45 L=173)   361    122        0 (1749/ 801-1804/ 915)  39.8 W41-W48 S28-S22
      9  855  647 (S13E19 L=109)   226     41        0 ( 827/ 629- 885/ 665)   0.0 E21-E17 S13-S11
     10 1351  571 (S31W08 L=136)   513     74        1 (1307/ 536-1389/ 611)   8.1 W04-W11 S32-S28
     11 1068  436 (S31E14 L=114)   306     60        0 (1051/ 404-1083/ 467)  68.2 E15-E12 S33-S28
     12  824  377 (S29E32 L= 96)   350    134        0 ( 810/ 321- 836/ 430)  45.0 E35-E28 S32-S25
  

    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