Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1112 1520 (N27E09 L= 33)   436     76        1 (1086/1485-1143/1551) -26.6 E11-E05 N24-N28
      2  665 1373 (N07E30 L= 12)   298     96        0 ( 612/1351- 704/1394)  45.0 E32-E28 N04-N08
      3 1846  985 (N17W57 L= 99)  2638    279        3 (1782/ 840-1879/1104) -42.9 W47-W62 N08-N22
      4 1647  984 (N10W39 L= 81)   355    120        1 (1608/ 928-1671/1042) -26.6 W34-W42 N07-N12
      5  724  971 (S15E16 L= 26)   391    125        1 ( 662/ 924- 780/1004)  54.5 E18-E13 S18-S11
      6  331  938 (S26E44 L=358)   715    152        2 ( 278/ 917- 386/ 976)  29.1 E48-E39 S27-S22
      7  726  804 (S24E12 L= 30)   408    110        1 ( 668/ 780- 781/ 832)   7.1 E16-E08 S25-S23
      8  418  699 (S39E34 L=  8)  1392    229        3 ( 331/ 612- 495/ 786) -25.0 E42-E26 S43-S34
      9 1316  716 (S14W25 L= 67)   539    154        1 (1283/ 637-1367/ 776)  90.0 W23-W27 S18-S08
     10  630  630 (S38E15 L= 27)   562    126        0 ( 568/ 612- 699/ 654)   5.2 E20-E09 S39-S35
     11 1107  628 (S24W16 L= 58)   360    129        0 (1040/ 613-1182/ 641)  18.4 W12-W21 S26-S22
     12 1534  396 (S25W55 L= 97)  1376    288        0 (1398/ 348-1646/ 446)  42.9 W48-W62 S31-S18
  

    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