Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1206 1831 (N52W52 L= 60)   453     45        2 (1179/1806-1231/1849)  -8.1 W47-W57 N50-N53
      2  302 1281 (N36E51 L=317)   865    181        2 ( 252/1221- 362/1322) -25.0 E58-E43 N31-N39
      3 1042 1290 (N22W07 L= 15)   295    134        0 ( 981/1252-1085/1315) -45.0 W03-W09 N18-N24
      4  928 1129 (N15E03 L=  5)  1025    254        3 ( 886/1037- 989/1253)  85.9 E04-E01 N08-N22
      5  438 1113 (N24E37 L=331)   491    147        0 ( 391/1072- 511/1151) -33.7 E40-E31 N20-N26
      6  626 1038 (N16E24 L=344)   270    133        1 ( 583/ 993- 672/1067)  18.4 E27-E20 N13-N17
      7  151  928 (N16E71 L=297)  1219     97        2 ( 100/ 902- 184/ 965)   8.5 E84-E64 N14-N17
      8 1639  887 (S16W36 L= 44)   451    116        1 (1615/ 834-1662/ 940)  71.6 W33-W38 S18-S12
      9  301  836 (N09E53 L=315)   905    153        3 ( 259/ 807- 358/ 898) -45.0 E57-E47 N07-N12
     10 1227  606 (S22W03 L= 11)   315     76        0 (1194/ 569-1256/ 642)  31.0 W00-W05 S24-S20
     11  844  343 (S32E31 L=337)   618     94        1 ( 824/ 296- 870/ 387)  45.0 E34-E27 S35-S28
     12  672  296 (S32E47 L=321)  1291    179        2 ( 634/ 219- 706/ 372)  57.3 E53-E43 S39-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