Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  976 1917 (N64W46 L= 37)   514     55        0 ( 950/1888- 992/1939)  23.2 W39-W53 N60-N66
      2  919 1892 (N64W35 L= 26)   452     45        1 ( 898/1864- 933/1913)  36.9 W30-W41 N60-N66
      3  788 1827 (N64W13 L=  4)  1685    253        1 ( 692/1747- 889/1894)   8.1 E02-W33 N59-N67
      4  922 1620 (N42W11 L=  2)  1408    310        0 ( 784/1592-1078/1651) -21.8 W00-W20 N35-N45
      5 1683 1544 (N15W61 L= 51)   471     55        1 (1648/1525-1708/1566)   0.0 W56-W64 N14-N16
      6 1218 1507 (N26W25 L= 16)   386     84        0 (1199/1463-1242/1543)  26.6 W22-W28 N23-N28
      7 1836 1427 (N05W70 L= 61)  1655    117        5 (1808/1373-1867/1478)  74.1 W66-W76 N00-N07
      8  383 1420 (N43E38 L=313)   508     60        0 ( 354/1396- 426/1437)  -7.1 E41-E33 N41-N44
      9 1509 1309 (N08W36 L= 26)   660    154        0 (1469/1234-1545/1385) -90.0 W34-W37 N02-N12
     10  159 1037 (N23E64 L=287)  4705    303        9 (  94/ 940- 225/1143)  80.5 E75-E55 N17-N29
     11 1184 1092 (N04W10 L=  1)   440    120        0 (1141/1040-1228/1143) -74.1 W08-W11 N00-N07
     12  267 1067 (N23E50 L=301)   370     63        1 ( 239/1035- 286/1095)  71.6 E52-E48 N21-N24
     13  406  914 (N12E40 L=310)   934    213        1 ( 333/ 861- 482/ 993) -51.3 E45-E35 N07-N17
     14 1180  795 (S12W03 L=354)   517    154        0 (1151/ 727-1210/ 877)  49.4 W00-W06 S15-S08
     15  142  745 (N06E75 L=276)  2269    126        1 ( 115/ 689- 175/ 808) -76.0 E82-E70 N02-N10
     16 1747  592 (S42W44 L= 35)   720    109        1 (1695/ 572-1797/ 616) -18.4 W37-W52 S45-S38
     17  564  534 (S13E41 L=310)   606    150        1 ( 534/ 465- 604/ 602)  28.6 E46-E35 S16-S09
     18 1579  513 (S42W26 L= 16)   500     83        0 (1543/ 483-1622/ 545)   8.1 W22-W29 S43-S40
  

    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