Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  819 1791 (N61W18 L= 18)   518     43        0 ( 791/1766- 844/1811)  78.7 W15-W21 N58-N63
      2 1224 1782 (N43W47 L= 47)   855    144        1 (1148/1753-1295/1803) -48.8 W43-W50 N39-N47
      3 1564 1748 (N28W73 L= 73)   761     79        1 (1520/1734-1607/1765) -18.4 W68-W80 N26-N30
      4  616 1659 (N58E11 L=348)  1272    225        3 ( 514/1594- 700/1717)   0.0 E24-W00 N55-N61
      5 1365 1611 (N28W42 L= 41)   770    180        0 (1286/1591-1460/1631) -32.5 W37-W48 N24-N31
      6 1574 1556 (N19W54 L= 53)   571    102        1 (1539/1528-1616/1580) -59.0 W50-W57 N15-N20
      7  430 1543 (N55E35 L=324)   478     59        1 ( 407/1502- 455/1572)  53.1 E38-E31 N52-N57
      8 1472 1423 (N15W39 L= 38)   813    293        0 (1333/1366-1588/1465) -45.0 W31-W44 N08-N21
      9  860 1352 (N29E01 L=359)   213     77        1 ( 833/1331- 896/1378) -63.4 E01-W01 N26-N31
     10  433 1174 (N29E35 L=325)   339    112        0 ( 378/1151- 482/1197)   0.0 E39-E30 N28-N30
     11  611 1055 (N18E23 L=336)   206     31        1 ( 593/1037- 636/1074) -18.4 E24-E21 N16-N19
     12 1629  711 (S30W30 L= 29)   680    159        0 (1558/ 687-1681/ 762)  -5.7 W24-W34 S31-S27
     13 1086  601 (S19E08 L=351)   426    153        0 (1071/ 533-1105/ 676)  66.0 E10-E04 S24-S15
     14  897  439 (S24E25 L=334)   763    284        0 ( 804/ 358- 945/ 558)  26.6 E34-E17 S28-S17
     15 1469  443 (S43W13 L= 12)  2004    343        4 (1326/ 373-1604/ 538)   2.3 W00-W26 S45-S39
     16 1053  402 (S31E17 L=342)   367    104        0 (1020/ 353-1085/ 439)  26.6 E21-E13 S34-S29
     17  641  300 (S27E51 L=308)   952    141        2 ( 579/ 266- 700/ 344) -63.4 E53-E48 S31-S22
     18 1101  255 (S45E23 L=336)  1510    221        1 (1034/ 214-1192/ 320) -82.4 E26-E15 S51-S36
     19  524  293 (S25E64 L=295)   914     78        2 ( 486/ 270- 558/ 313) -31.0 E68-E58 S27-S23
  

    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