Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1621 1591 (N22W55 L=255)  1513    261        2 (1523/1552-1719/1649) -28.1 W49-W64 N18-N28
      2 1710 1436 (N11W56 L=255)   701    165        0 (1638/1373-1753/1473) -45.0 W50-W58 N07-N15
      3 1161 1423 (N22W17 L=216)  2330    335        6 ( 924/1387-1443/1446) -16.9 W01-W34 N16-N27
      4 1263 1261 (N11W19 L=218)   345    139        0 (1196/1226-1312/1284) -29.7 W14-W21 N07-N12
      5 1579 1256 (N04W39 L=238)   270     34        0 (1558/1240-1606/1275) -45.0 W37-W40 N02-N05
      6  145 1211 (N29E62 L=138)   744    218        0 ( 110/1189- 188/1248) -21.8 E66-E56 N27-N31
      7  720 1156 (N15E15 L=184)   790    198        2 ( 634/1070- 791/1216) -57.5 E20-E12 N09-N20
      8 1807 1050 (S12W53 L=252)  1972    350        9 (1743/ 966-1869/1124)  26.6 W46-W61 S16-S08
      9  217 1080 (N20E54 L=145)   562     88        2 ( 184/1043- 245/1106) -45.0 E57-E51 N17-N21
     10 1547  951 (S13W30 L=230)   279    113        1 (1499/ 908-1567/ 978) -68.2 W27-W31 S15-S10
     11 1836  903 (S22W53 L=253)   680    131        1 (1797/ 841-1859/ 966)  45.0 W49-W56 S24-S18
     12  296  900 (N08E50 L=149)   505    208        0 ( 231/ 884- 370/ 913)  -8.7 E56-E43 N06-N09
     13  255  617 (S08E64 L=135)   611    108        1 ( 238/ 569- 272/ 671)  71.6 E66-E62 S10-S04
     14  306  610 (S09E59 L=141)   524    113        0 ( 292/ 564- 318/ 657)  39.8 E61-E55 S11-S06
     15 1214  508 (S33W01 L=201)   919    213        0 (1107/ 445-1277/ 595)  24.8 E07-W06 S35-S28
     16  642  506 (S21E36 L=164)   560    177        1 ( 591/ 452- 674/ 580)  31.0 E41-E31 S23-S17
     17  755  407 (S30E32 L=167)  1141    243        2 ( 697/ 328- 817/ 496)  29.4 E40-E24 S34-S25
     18  292  456 (S18E70 L=130)   925     82        0 ( 256/ 418- 321/ 490) -51.3 E72-E67 S20-S15
  

    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