Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  956 1801 (N39E35 L=348)   666     90        2 ( 919/1766- 992/1824) -14.0 E38-E30 N37-N41
      2  487 1649 (N16E57 L=325)  1498    250        1 ( 450/1548- 544/1753) -56.3 E64-E51 N09-N24
      3  586 1502 (N09E41 L=341)  4314    703        4 ( 488/1371- 761/1714) -61.4 E54-E31 N00-N22
      4 1374 1558 (N36W05 L= 27)   454     89        1 (1353/1519-1393/1605) -45.0 W02-W07 N33-N39
      5 1392 1478 (N31W09 L= 31)   441     48        0 (1362/1457-1425/1499)  11.3 W06-W11 N29-N32
      6 1047 1396 (N15E09 L= 13)   339    136        0 (1004/1336-1104/1437) -20.6 E12-E04 N13-N16
      7  213 1274 (S12E65 L=317)  1716    186        1 ( 127/1242- 269/1322)  10.3 E81-E59 S13-S07
      8  201 1156 (S19E63 L=319)   420     45        1 ( 185/1129- 215/1182)  76.0 E65-E60 S20-S16
      9 1340 1137 (N09W14 L= 36)   304     53        1 (1324/1102-1354/1175) -78.7 W12-W15 N06-N11
     10  437 1121 (S16E39 L=344)   307     68        0 ( 406/1099- 461/1150)  53.1 E40-E37 S17-S13
     11  187  975 (S30E62 L=320)   909    110        1 ( 135/ 928- 229/1004)  21.8 E71-E56 S33-S27
     12  169  923 (S33E65 L=317)   739    103        0 ( 131/ 883- 214/ 950)  21.8 E74-E58 S36-S30
     13  725  767 (S29E10 L= 12)   664    150        2 ( 651/ 741- 794/ 799)  33.7 E15-E06 S32-S25
     14 1420  702 (S12W33 L= 55)   430    184        0 (1335/ 680-1506/ 730)  18.4 W26-W38 S13-S08
     15 1232  570 (S25W27 L= 49)   456    164        1 (1156/ 558-1319/ 589)  26.6 W22-W32 S26-S21
     16  971  539 (S35W13 L= 35)   271    102        0 ( 921/ 524-1015/ 561)  15.9 W08-W15 S36-S33
     17 1510  518 (S18W48 L= 70)  1048    158        0 (1441/ 492-1581/ 557)  41.2 W42-W51 S21-S14
     18 1123  403 (S38W30 L= 52)  1930    344        1 ( 994/ 327-1272/ 480)  69.8 W24-W34 S47-S28
     19  943  379 (S46W20 L= 42)   264     36        2 ( 922/ 355- 957/ 405) -53.1 W17-W22 S48-S44
  

    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