Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1167 1819 (N52W23 L= 63)   484    100        1 (1146/1786-1196/1846)  40.6 W19-W26 N48-N54
      2 1632 1421 (N19W45 L= 86)   856    186        4 (1576/1367-1687/1465) -35.5 W40-W48 N15-N21
      3  759 1371 (N21E14 L= 26)   363    106        1 ( 715/1351- 813/1403) -23.2 E17-E10 N19-N23
      4 1293 1299 (N13W19 L= 59)   245     69        0 (1261/1282-1335/1313)   0.0 W16-W21 N12-N13
      5  352 1077 (N06E44 L=357)   498    112        1 ( 337/1017- 373/1141) -90.0 E44-E42 N02-N10
      6  259 1027 (N04E52 L=349)  1407    257        1 ( 207/ 903- 295/1131) -69.0 E57-E49 S03-N10
      7  126  969 (N02E69 L=332)  5148    578       14 (  75/ 857- 167/1101) -77.9 E79-E63 S04-N10
      8 1747 1020 (S06W48 L= 89)   268     40        0 (1730/ 992-1765/1047)  45.0 W46-W50 S07-S04
      9  852  806 (S14E12 L= 29)   188     48        1 ( 839/ 772- 871/ 831) -90.0 E12-E10 S15-S11
     10  274  596 (S22E60 L=341)   510    150        4 ( 257/ 554- 292/ 650)  80.5 E62-E57 S24-S18
     11  798  366 (S42E25 L= 16)   820    160        2 ( 723/ 340- 871/ 409)   8.1 E31-E17 S44-S39
  

    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