Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1220 1858 (N51W50 L=117)   448     14        1 (1198/1839-1241/1875)  14.0 W45-W53 N49-N52
      2 1467 1487 (N19W39 L=107)  1596    377        2 (1299/1436-1605/1584)   0.0 W27-W55 N16-N22
      3  958 1520 (N33W07 L= 75)  1844    385        2 ( 807/1444-1132/1574) -35.8 E02-W16 N24-N38
      4  342 1413 (N40E45 L= 23)  2587    300        7 ( 257/1307- 406/1556) -35.0 E69-E36 N32-N53
      5  369 1085 (N18E41 L= 27)  1481    238        3 ( 316/ 990- 423/1184) -67.4 E45-E37 N12-N24
      6  118 1050 (N21E72 L=356)  3453    244        8 (  90/ 945- 151/1155)  90.0 E79-E66 N14-N27
      7 1282  848 (S13W11 L= 79)   416    123        0 (1253/ 786-1306/ 916)  60.3 W08-W13 S16-S09
      8 1745  799 (S27W45 L=113)   396     51        1 (1718/ 774-1776/ 815) -18.4 W42-W48 S28-S25
      9  990  688 (S16E09 L= 59)   509    153        0 ( 934/ 638-1057/ 737) -58.0 E11-E06 S20-S12
     10  449  493 (S17E52 L= 16)   574    110        0 ( 423/ 465- 493/ 534) -78.7 E54-E48 S19-S14
     11  631  410 (S26E42 L= 26)   430    107        2 ( 601/ 369- 652/ 463)  40.6 E46-E39 S28-S22
     12  441  400 (S23E60 L=  8)   526     85        1 ( 419/ 358- 459/ 440)  50.2 E64-E57 S26-S20
  

    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